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(54) Title: METHOD FOR COMPUTING FAST FOURIER TRANSFORM AND INVERSE FAST FOURIER TRANSFORM 

(54) Titre: PROCEDE DE CALCUL DE LA TRANSFORMEE DE FOURIER RAPIDE ET DE LA TRANSFORMEE DE FOURIER 
RAPIDE INVERSE 

(57) Abstract 

The invention concerns 
a method for computing the 
fast Fourier transform and the 
inverse fast Fourier transform 
of a series of N real sam- 
ples x(n), with N power of 2, 
functioning on the basis of an 
algorithm with time interlace 
and delivering the series of 
samples X(n) in ascending or- 
der of the index n and which 
uses reduced computing and 
storing means. The inven- 
tion also concerns a method 
for computing the fast Fourier 
transform and the inverse fast 
Fourier transform of a se- 
ries of N conjugated complex 
samples X(n), with N power 
of 2, functioning on the basis 
of an algorithm with frequen- 
tial interlace. The invention 
is useful for treating images or acoustic 
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signals and for multicarrier modulation. 



(57) Abrege* 

L* invention conceme un proc^de* de calcul de la transformer de Fourier rapide ou de la transformee de Fourier rapide inverse d*une 
suite de N echantillons reels x(n), avec N puissance de 2, fonctionnant selon un algorithme avec entrelacement temporel et dflivrant la 
suite des echantillons X(n) dans l'ordre croissant de Tindice n et qui utilise des moyens de calcul et de stockage nfduits. L* invention 
concerne egalement un proc&16 de calcul de la transformee de Fourier rapide ou de la transformee de Fourier rapide inverse d'une suite 
de N echantillons complexes conjugues X(n), avec N puissance de 2, fonctionnant selon un algorithme avec entrelacement frequentiel. 
Application: traitement d 'images ou de signaux acoustiques, modulation multiporteuse. 
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PROCEDE DE CALCUL DE LA TRANS FORME E DE FOURIER RAP IDE 
ET DE LA TRANSFORMEE DE FOURIER RAPIDE INVERSE 

L 1 invention concerne un procede die calcul de la 

5 transformee de Fourier rapide ou de la transformee de 

Fourier rapide inverse d'une suite de nombres reels ou 
d'une suite d 1 echantillons complexes conjugues. 

La transformee de Fourier est sans doute l'un des 
D outils les plus importants dans 1' analyse, la 
conception et la mise en oeuvre des algorithmes 
relevant du traitement du signal, et 1' existence 
d 1 algorithmes efficaces tels que celui de la 
transformee de Fourier rapide a largement contribue & 
5 cette situation. Quoique la plupart des algorithmes de 
transformee de Fourier soient con?us pour transformer 
des suites de nombres complexes, il existe cependant de 
nombreuses applications, telles que le traitement 
d 1 images ou de signaux acoustiques ou certains types de 
0 la modulation multiporteuse, dans lesquelles les suites 
a transformer sont des nombres reels, 

De fa<?on generale, la transformee de Fourier directe et 
la transformee de Fourier inverse Stablissent 
5 respectivement entre deux suites de N nombres 
complexes, x(n) et X(n) , les relations suivantes : 



kn 



avec //€[0...AT-l] et w kn = e 




*(//) = — £ X(k)w' kn avec ?i e[0...N-\] 
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En 19 65, J.W Cooley et J.W Tukey ont deer it un 
algorithme permettant de calculer rapideitient la 
transf ormee de Fourier d'une suite de N nombres 
complexes dans un article intitule "An Algorithm for 
5 the Machine Calculation of Complex Fourier Series", 
Math. Computation, Vol. 19, 1965, pp. 297-3 01. Cet 
algorithme est tres interessant lorsque N est une 
puissance de deux car il est alors particulierement 
simple a mettre en oeuvre. Cet algorithme requiert /z 

10 etapes de calcul ou /z=log 2 (N) . II repose sur une 
decomposition de la suite a transformer en deux sous- 
suites entrelacees. II existe deux types 
d 1 entrelacement : 1 1 entrelacement temporel et 
1 1 entrelacement frequentiel. Ces deux types 

15 d 1 entrelacement seront exposes plus en detail dans la 
suite de la description en reference aux figures 1A et 
IB. 

Des algorithmes particuliers ont ete developpes a 
20 partir de cet algorithme de base pour traiter le cas 
des suites de nombres reels. La transformation de 
Fourier d'une suite de 8 nombres reels selon un 
algorithme avec entrelacement temporel et selon un 
algorithme avec entrelacement frequentiel est illustree 
2 5 aux figures 1A et IB. Pour chaque algorithme de 
transf ormee de Fourier avec entrelacement temporel, il 
existe un algorithme avec entrelacement frequentiel qui 
correspond a une double inversion de la suite des 
operations de transformation d'une part, et, pour 
30 chaque circuit croisillon, des operations de 
transformation en elles-memes d ' autre part. Quel que 
soit 1 1 entrelacement choisi, le procede de 
transformation necessite trois etapes de transformation 
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EO, El et E2, lesquelles etapes sont mises en oeuvre 
par un jeu de quatre circuits croisillons CC, 
communement appeles "butterfly" en langage technique. 
Chaque circuit croisillon, represents sur la figure par 
5 un point de croiseraent entre deux colonnes de nombres, 
effectue des operations de calcul sur deux nombres, 
reels ou complexes* Les symboles R et C identifient 
respectivement un nombre reel et un nombre complexe. 
L ' ordre des nombres reels et/ou complexes a l 1 issue des 
10 etapes de transformation est dependant de 
1 1 entrelacement choisi . 

L' algorithme avec entrelacement temporel est 
generalement choisi pour calculer la transformee de 

15 Fourier d f une suite de nombres reels en raison de la 
repartition symetrique des nombres reels et complexes a 
travers les etapes. En revanche, l'algorithme avec 
entrelacement frequentiel est plus adapte pour la 
transformation directe ou inverse d'une suite nombres 

2 0 complexes conjugues. 

Lorsque la suite a transformer x(n) est reelle, la 
transformee de Fourier verifie la relation suivante: 

x(n) est reelle si et seulement si 
25 X(n) = X*(-n)= X*(N-n) ; (1) 
ou * designe 1' operation de conjugaison. 

Pour une suite x(n) de N nombres reels, on deduit de 
cette relation les resultats suivants: 
- X(0) et X(N/2) sont reels; 
30 - X(n) = X*(N-n) pour 1 < n < N/2 - 1 

La relation (1) met en evidence la presence 
d' informations redondantes dans la suite X(n) . 
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II convient de noter que le procede de transformation 
est generalement mis en oeuvre par un unique jeu de 
circuits croisillons dont le mode operatoire est 
modifie au fur et a mesure de la transformation. A 
5 chaque changement de mode operatoire, les resultats 
sont stockes dans une memoire comportant N emplacements 
memoire, les echantillons de sortie d'un circuit 
croisillon se substituant dans la memoire aux 
echantillons d 1 entree de meme rang correspondants . 

10 Cette methode d 1 application de 1'algorithme est 
couramment appelee methode "in place". Cette methode 
presente un avantage majeur : si les elements de la 
suite x(n) sont traites a la premiere etape de 
transformation dans l f ordre binaire inverse de 1 ' indice 

15 n ("bit-reversed"), les nombres de la suite X(n) sont 
delivres a la derniere etape de transformation dans 
I'ordre croissant de !• indice n et vice versa. 

Un procede de transformation connu est presente a titre 
20 d'exemple a la figure 2. Ce procede effectue la 
transformation de Fourier d'une suite reelle x(n) selon 
un algorithme complexe avec entrelacement temporel. 
Dans cet exemple, la suite x(n) a transformer comporte 
seize echantillons reels, x(0) a x(15). Le procede de 
25 transformation comporte quatre etapes de transformation 
Ep avec 0<p<3. Les echantillons de la suite x(n) sont 
presentes a la premiere etape de transformation dans 
I'ordre binaire inverse de leur indice n 
("bit-reversed") . 

30 

A ce stade des explications, il convient de definir les 
termes employes dans la suite de la description. Le 
rang d'un echantillon s'entend de la position occupee 
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par celui-ci dans la suite des echantillons a laquelle 
il appartient. L'indice d'un echantillon correspond 
alors au rang de depart de cet echantillon. 



5 Les resultats intermediaires obtenus aux differentes 
etapes de transformation sont representes par les 
suites A(n) , B(n) et C(n) . Les echantillons des suites 
x(n), A(n) , B(n) , C(n) et X(n) sont stockes dans des 
doubles emplacements memoire, un emplacement memoire 

10 etant reserve pour la partie reelle de 1 1 echantillon et 
1' autre emplacement etant reserve pour sa partie 
imaginaire. A R (n) et Aj(n) designent respectivement la 
partie reelle et la partie imaginaire de 1 1 echantillon 
d'indice n de la suite A(n) . Les circuits croisillons 

15 sont representes sur la figure par des points de 
croisement entre des colonnes d 1 emplacements memoire. A 
chaque circuit croisillon est affecte un coefficient W s 
symbolise sur la figure 2 par un couple de coordonnees 
A/B ou A et B designent respectivement la partie reelle 

20 et la partie imaginaire du coefficient W s . Les 
coordonnees 1/0 et 0/-1 sont attributes respectivement 
aux coefficients W°=l et W N / 4 =W 4 =-j. Par souci de 
clarte et pour simplifier leur formulation, les autres 
coefficients W s ont ete representes par les couples 

25 suivants: 

W 1 -> 2 / -4 W 5 -> -4 / -2 

W 2 -> 3 / -3 W 6 -> -3 / -3 

W 3 -> 4 / -2 W 7 -> -2 / -4 



3 0 Ces couples de coordonnees sont representes 

graphiquement a la figure 3. Les coordonnees A et B 

representent en realite respectivement une valeur de 

cosinus et une valeur de sinus. Ce coefficient W s 
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intervient dans le calcul effectue par le circuit 
croisillon. Par ailleurs, les circuits croisillons sont 
repartis a chaque etape de transformation dans N/2P +1 
blocs de calcul, chaque bloc de calcul comportant 2^ 
5 circuits croisillons. Dans la suite de la description, 
le parametre q designe le rang des blocs de calcul dans 
une meme etape de transformation; q est compris entre 0 
et (N/2P +1 )-1. 

A la premiere etape de transformation E Q , les circuits 
croisillons sont repartis en huit blocs de calcul 
comprenant chacun un circuit croisillon effectuant une 
operation sur deux echantillons complexes ou reels. Si 
on designe par el et e2 les echantillons appliques sur 
les entrees d'un circuit croisillon, ce dernier delivre 
en sortie des echantillons si et S2 definis de la 
maniere suivante : 

si = el + W s .e2 et s2 = el -(W s .e2) 

ou W s est le coefficient affecte audit croisillon. 

Pour cette premiere etape de transformation, le 
coefficient W°=l est affecte aux huit circuits 
croisillons. Etant donne que les echantillons x(n) et 
le coefficient W° sont reels, les echantillons A(n) 
obtenus en fin d' etape E Q sont reels. 

Pour la seconde etape de transformation, E lf les 
circuits croisillons sont repartis en quatre blocs de 
calcul comprenant chacun deux circuits croisillons. Le 
30 coefficient W°=l est affecte au premier de ces circuits 
croisillons; ainsi, le premier circuit croisillon de 
chaque bloc de calcul delivre deux echantillons reels. 
Le deuxieme circuit croisillon des blocs de calcul est 
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associe au coefficient W N ^ 4 =W 4 =-j et genere deux 
echantillons complexes conjugues. Les echantillons de 
sortie obtenus en fin d* etape E^^ sont designes par la 
suite B (n) . 

5 

Pour la troisieme etape de transformation, E 2 , les 
circuits croisillons sont repartis dans deux blocs de 
calcul comprenant chacun quatre circuits croisillons 
associes respectivement aux coefficients W°, W 2 ' W 4 et 

10 W 6 . Les echantillons de sortie a 1* etape E 2 sont 
designes par la suite C(n) . Enf in, pour la quatrieme 
etape de transformation, E 3 , il est prevu un unique 
bloc de calcul comportant huit circuits croisillons 
associes respectivement aux coefficients W° , W 1 , W 2 , 

15 W 3 , W 4 , W 5 , W 6 et W 7 . Cette etape de transformation 
genere la suite transformee X(n) . 

Etant donne la relation (1) , la suite X(n) comporte 
d'une part des echantillons reels, X(0) et X(8) , et 

20 d* autre part des echantillons complexes, X(l) a X(7) et 
X(9) a X(15), les echantillons X(15) a X(9) etant 
respectivement les conjugues des echantillons X(l) a 
X(7) . La suite X(n) contient done des informations 
redondantes. Les emplacements memoire dessines en trait 

25 epais sur la figure 2 designent les emplacements 
memoire renfermant les valeurs conjuguees des 
echantillons complexes contenus dans les emplacements 
memoire qui leur sont associes par une fleche. Les 
suites de resultats intermediaires B(n) et C(n) 

3 0 contiennent egalement des informations redondantes. 

II est alors possible de supprimer ces informations 
redondantes afin de reduire de moitie la taille de la 
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memoire de stockage des echantillons ainsi que le 
nombre des circuits croisillons. 

Cependant, une suppression des informations redondantes 
5 stockees dans les emplacements memoire dessines en 
trait gras sur la figure 2 implique de reorganiser 
totalement les etapes de transformation de la figure 2. 
Cette reorganisation de la transformation a pour effet 
de modifier 1 1 ordre de sortie des echantillons X(n) . 

10 

Le probleme est done de reduire la taille de la memoire 
de stockage et le nombre des circuits croisillons tout 
en conservant 1' ordre de sortie des echantillons X(n) . 
L* invention a pour but de proposer un procede de calcul 

15 de la transformee de Fourier rapide ou de transformee 
de Fourier rapide inverse d'une suite de N echantillons 
reels x(n) , avec N puissance de 2, fonctionnant selon 
un algorithme avec entrelacement temporel, qui delivre 
la suite des echantillons X(n) dans 1 1 ordre croissant 

20 de I'indice n et utilise des moyens de calcul et de 
stockage reduits. 

Dans ce but, l 1 invention a pour objet un procede de 
calcul de la transformee de Fourier rapide ou de la 

2 5 transformee de Fourier rapide inverse d'un signal 
numerique defini par une suite de N echantillons reels 
de depart x(n), avec N puissance de deux et n e 
[0..N-1], comportant des etapes successives de 
transformation pour transformer des echantillons 

30 d 1 entree en echantillons de sortie, 1» ensemble des 
etapes de transformation etant effectue au moyen d'un 
unique jeu de circuits croisillons a plusieurs entrees 
et plusieurs sorties dont le mode operatoire est change 
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selectivement a chaque etape de transformation, les 
echantillons d 1 entree et de sortie de chaque etape de 
transformation etant stockes dans une memoire de 
stockage, une suite de N echantillons de sortie y(n) 
5 representatif s de la transformee de Fourier rapide ou 
rapide inverse des echantillons de depart x(n) etant 
delivree a la derniere etape de transformation, 

caracterise en ce que les echantillons de sortie 
y(n) sont reels, 

10 et en ce que les echantillons de sortie d*un 

circuit croisillon se substituent dans la memoire de 
stockage aux echantillons d 1 entree de meme rang 
correspondants, de telle sorte que, si les echantillons 
de depart x(n) traites a la premiere etape de 

15 transformation sont classes dans 1 1 ordre binaire 
inverse de leur indice n, les echantillons de sortie 
y(n) sont delivres a la derniere etape de 
transformation dans 1 1 ordre croissant de 1 • indice n, 
lesquels echantillons de sortie sont definis par les 

20 relations suivantes: 
y(0) = Re[X(0)] 

y(n) = Re[X( (n+1) /2) ] pour n impair et 

different de N-l 
y(n) = Im[X(n/2)] pour n pair et different 

25 de 0 

y(N-l) = Re[X(N/2)] 
ou les echantillons X(n), avec n e [0..N-1], designent 
les echantillons complexes de la suite correspondant a 
la transformee de Fourier rapide ou rapide inverse de 
30 la suite des echantillons de depart x(n) . 

Pour les precedes fonctionnant selon un algorithme avec 
entrelacement f requentiel , 1 1 invention concerne 
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egalement un procede de calcul de la transformee de 
Fourier rapide ou de la transformee de Fourier rapide 
inverse d'un signal numerique defini par une suite de N 
echantillons complexes X(n) conjugues deux a deux 
representes par une suite de N echantillons reels de 
depart y(n) , avec N puissance de deux et n e [0..N-1], 
les echantillons de depart y(n) etant definis de la 
maniere suivante : 

y(0) = Re[X(0) ] 

y(n) = Re[X((n+l)/2)] pour n impair et 

different de N-l 

y(n) = Im[X(n/2) ] pour n pair et different 

de 0 

y(N-l) = Re[X(N/2) ] 
lequel procede de calcul comporte des etapes 
successives de transformation pour transformer des 
echantillons d ' entree en echantillons de sortie, une 
suite de N echantillons de sortie x(n) representatif s 
de cette transformee de Fourier rapide ou rapide 
inverse etant delivree a la derniere etape de 
transformation, l 1 ensemble des etapes de transformation 
etant effectue au moyen d'un unique jeu de circuits 
croisillons a plusieurs entrees et plusieurs sorties 
dont le mode operatoire est change selectivement a 
chaque etape de transformation, les echantillons 
d' entree et de sortie de chaque etape de transformation 
etant stockes dans une memoire de stockage, 

caracterise en ce que les echantillons de sortie 
x(n) sont reels, 

et en ce que les echantillons de sortie d'un 
circuit croisillon se substituent dans la memoire de 
stockage aux echantillons d' entree de meme rang 
correspondants, de telle sorte que, si les echantillons 
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de depart y(n) traites a la premiere etape de 
transformation sont classes dans 1 1 ordre croissant de 
l'indice n, les echantillons de sortie x(n) sont 
delivres a la derniere etape de transformation dans 
5 1 1 ordre binaire inverse de l'indice n. 

Les procedes de calcul de 1 ' invention executent des 
operations sur des echantillons reels et utilisent en 
consequence des moyens de calcul et de stockage reduits 
10 par rapport au precede presente a la figure 2. 

Selon une autre caracteristique de l 1 invention, les 
circuits croisillons transf orment , a chaque etape de 
transformation, des couples d 1 echantillons d* entree, 

15 les rangs des echantillons d 1 entree d'un meme couple 
dans la suite des echantillons d 1 entree de ladite etape 
de transformation etant symetriques par rapport a un 
milieu entre les valeurs de rang extremes des 
echantillons d' entree transf ormes par ledit circuit 

20 croisillon. Les echantillons d" entree traites par un 
meme circuit croisillon sont ainsi relies deux a deux 
par une symetrie. II en resulte une gestion simplifiee 
de l'adressage des echantillons. 

2 5 Selon un autre aspect de l 1 invention, le procede 
comporte de preference jz-l etapes de transformation E n 

XT 

avec ju=log 2 (N) et p € [0../X-2]. 

D'autres caracteristiques et avantages de 1* invention 
30 apparaitront a la lecture de la description detaillee 
qui suit et qui est faite en reference aux dessins 
annexes dans lesquels : 
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les figures 1A et IB, deja decrites, 
representent respectivement une transformation de 
fourier de huit nombres reels selon un algorithme avec 
entrelacement temporel et selon un algorithme avec 
5 entrelacement frequentiel; 

la figure 2, deja decrite, illustre la 
transformation d'une suite de 16 nombres reels en une 
suite de 16 nombres complexes selon un algorithme 
complexe avec entrelacement temporel; 
1° - ia figure 3, deja decrite, represente 

graphiquement la correspondance entre les coefficients 
W s et des couples de coordonnees A/B; 

la figure 4 represente des modifications 
appliquees sur une partie de la transformation de la 
15 figure 2; 

- la figure 5 illustre une transformation modifiee 
traitant uniquement des nombres reels; 

- la figure 6 est une representation du procede de 
calcul de la transformee de Fourier rapide selon 

2 0 1' invention; 

- les figures 7 A et 7B illustrent respectivement 
des permutations effectuees sur les circuits 
croisillons de bord et sur les circuits croisillons 
internes de rang impair de la transformation de la 

25 figure 5; 

- la figure 8 represente un mode de realisation du 
procede de transformation selon 1' invention, comportant 
/x etapes de transformation; 

- la figure 9 illustre un regroupement de circuits 
30 croisillons de bord selon un premier mode de 

realisation d'un procede de transformation comprenant 
M"l etapes de transformation; 
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- la figure 10 represente un premier mode de 
realisation d'un procede de transformation comprenant 
li-l etapes de transformation; 

- la figure 11 represente un modele de circuit 
5 croisillon se rapportant au procede de transformation 

illustre figure 10; 

- la figure 12 represente une variante du mode de 
realisation de la figure 10; 

- la figure 13 represente un modele de circuit 
10 croisillon se rapportant au mode de realisation de la 

figure 12; 

- la figure 14 illustre un regroupement de 
circuits croisillons de bord selon un deuxieme mode de 
realisation d'un procede de transformation comprenant 

15 p-1 etapes de transformation; 

- la figure 15 represente un deuxieme mode de 
realisation d'un procede de transformation comprenant 

etapes de transformation avec \i pair; 

- la figure 16 represente une variante du mode de 
20 realisation precedent avec fi impair; 

- la figure 17 represente un modele de circuit 
croisillon se rapportant aux modes de realisation 
presentes aux figures 15 et 16; 

- la figure 18 represente un troisieme mode de 
25 realisation d'un procede de transformation comprenant 

Ax-1 etapes de transformation; 

- la figure 19 represente un modele de circuit 
croisillon se rapportant au mode de realisation de la 
figure 18; 

3 0 - la figure 2 0 represente les adresses qui sont 

associees aux differents circuits croisillons mis en 
oeuvre dans le mode de realisation presente a la figure 
12; 



WO 99/30251 14 PCT/FR98/02636 



- la figure 21 represente les adresses qui sont 
associees a une partie des echantillons d*un procede de 
transformation traitant une suite de 3 2 echantillons 
reels; 

5 - la figure 22 represente un exemple de mode de 

realisation d'un procede de transformation fonctionnant 
selon un algorithme avec entrelacement frequentiel; 

- la figure 2 3 represente un modele de circuit 
croisillon se rapportant au mode de realisation de la 

10 figure 22 ; 

- la figure 24 represente les adresses qui sont 
associees aux differents circuits croisillons du mode 
de realisation illustre a la figure 22. 

15 Selon !• invention, on ne calcule qu'une partie des 
echantillons X(n), l 1 autre partie des echantillons 
etant redondante. On pourrait par exemple se limiter au 
calcul des nombres X(n) pour 0<n<N/2. Cependant, une 
solution plus adaptee consiste a calculer les nombres 

20 X(n) pour 0<n<N/4 et N/2<n<3N/4. Cette derniere 
solution est preferable car elle ne fait intervenir que 
les (N/4)+l premiers circuits croisillons pour le 
calcul des X(n) , les (N/4)-l derniers circuits 
croisillons pouvant alors etre supprimes. Les 

25 emplacements memoire ainsi liberes peuvent etre 
utilises pour memoriser la partie reelle ou la partie 
imaginaire des nombres X(n) restants. La taille de la 
memoire de stockage peut ainsi etre divisee par deux si 
on limite la taille des emplacements memoire au 

30 stockage d ! un nombre reel au lieu d'un nombre complexe. 

La partie reelle du nombre complexe X(n) est stockee 
dans 1 1 emplacement memoire qui lui a ete initialement 
affecte tandis que sa partie imaginaire est stockee 
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dans 1 •emplacement memoire affecte initialement au 
nombre X(N-n) • La meme operation peut etre effectuee 
pour les suites de resultats intermediates B(n) et 
C(n) . Le procede de transformation de l'invention se 
5 limite a calculer N echantillons reels y(n) au lieu de 
N echantillons complexes. Les echantillons reels y(n) 
sont definis de la maniere suivante: 
y(0) = Re[X(0) ] 

y(n) = Re[X( (n+1) /2) ] pour n impair et 
10 different de N-l 

y(n) = Im[X(n/2)] pour n pair et different 

de 0 

y(N-l) = Re[X(N/2)]. 

15 La suppression des informations redondantes et la 
reorganisation des emplacements memoire modifient 
fortement le schema de la figure 2. La figure 4 
illustre les reorganisations effectuees sur la partie 
se rapportant au calcul des echantillons A(l), A(3) , 

20 A<5), A(7), B(l), B(3), B(5), B(7), C(l), C(3), C(5) et 

C(7) . Les echantillons redondants a supprimer sont 
B(3), B(7), C(5) et C(7). Apres reduction de la taille 
des emplacements memoire et rearrangement des 
echantillons restants dans ces emplacements memoire, 

25 les circuits croisillons effectuent des calculs sur des 
nombres reels, Les circuits croisillons auxquels est 
affecte le couple 0/-1 effectuent des calculs sur deux 
nombres reels. En pratique, ils recopient sur leur 
premiere sortie le nombre present sur leur premiere 

3 0 entree et multiplient par -1 le nombre present sur leur 
seconde entree et le delivrent sur leur seconde sortie. 
Les circuits croisillons auxquels est affecte le couple 
1/0 effectuent une addition et une soustraction sur 
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deux nombres reels. Enfin, les autres circuits 
croisillons effectuent des operations sur quatre 
nombres reels* 

5 La transformation ainsi rearrangee est illustree a la 
figure 5. Sur cette figure, les circuits croisillons 
associes aux couples 1/0 et 0/-1 correspondant aux 
coefficients W° et W 4 sont connectes aux emplacements 
memoire par des traits epais. Cette figure montre que 

10 la reorganisation des etapes du procede de 
transformation modifient 1 1 ordre de sortie des 
echantillons y(n) et done 1 ' ordre de sortie de la suite 
X(n) . Par ailleurs, ce procede de transformation ne 
presente plus de symetrie particuliere permettant de 

15 relier deux a deux les echantillons traites par un meme 
circuit croisillon. II en resulte une gestion tres 
compliquee de 1'adressage des echantillons a appliquer 
sur les entrees des circuits croisillons, 

2 0 La figure 6 presente les etapes d'un procede de calcul 
de la transformee de Fourier rapide ou de la 
transformee de Fourier rapide inverse d'une suite de N 
nombres reels x(n) , avec N puissance de 2, fonctionnant 
selon un algorithme avec entrelacement temporel. II 

25 comporte essentiellement des etapes de transformation 2 
visant a transformer les N echantillons de depart x(n) 
classes dans 1 ' ordre binaire inverse de leur indice n 
en echantillons de sortie y(n) reels representatif s de 
cette transformee de Fourier classes dans 1' ordre 

30 croissant de 1' indice n. 

Avantageusement , il comporte une etape 
preliminaire 1 pour ranger les N echantillons reels de 
depart x(n) a transformer dans 1 ' ordre binaire inverse 



WO 99/30251 17 PCT/FR98/02636 



de leur indice n si les echantillons x(n) ne sont pas 
deja dans cet ordre et une etape finale 3 pour generer 
les N echantillons complexes X(n) correspondant a la 
transformee de Fourier rapide des echantillons de 
5 depart x(n) a partir des N echantillons reels y(n) 
obtenus a 1 9 issue des etapes de transformation. 

Les procedes qui seront detailles dans la suite de la 
description, seront plus particulierement destines a 

10 calculer la transformee de Fourier rapide d'une suite 
reelle. Aussi, les coefficients W s affectes aux 
circuits croisillons pour la mise en oeuvre du precede 
de l'invention seront du type e~3 ( 27rs / N ) avec s>0. Pour 
le calcul de la transformee de Fourier rapide inverse, 

15 le procede de calcul est identique, cependant le 
coefficient est du type (27TS/N) avec s >q. 

Pour obtenir simultanement des echantillons y(n) 
ordonnes dans 1' ordre croissant de 1 1 indice n et une 

2 0 symetrie de calcul, on propose selon 1' invention de 

modifier les operations de calcul effectuees par les 
circuits croisillons des blocs de calcul de rang q 
impair de la transformation illustree figure 5 
conformement aux schemas des figures 7 A et 7B. 

25 

S'agissant des circuits croisillons associes au 
coefficient 1/0 (figure 7A) des blocs de calcul de rang 
impair, on prevoit de permuter les deux sorties du 
circuit croisillon et de multiplier par -1 le resultat 

3 0 delivre sur la deuxieme sortie du circuit croisillon. 
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S'agissant des circuits croisillons a quatre entrees 
(figure 7B) , on prevoit de permuter les deux premieres 
sorties avec les deux dernieres. 

5 On applique cette methode a 1* ensemble de la 
transformation et on obtient alors un procede de 
transformation delivrant en sortie des echantillons 
y(n) dans I'ordre croissant de l'indice n. Ce procede 
est illustre a la figure 8. Le symbole ^) place au- 

10 dessus des blocs de calcul, designe les blocs de calcul 
dans lesquels ont ete modifies les circuits 
croisillons, c'est-a-dire les blocs de calcul de rang 
impair. Etant donne la symetrie de calcul du procede de 
transformation, les points de croisement representant 

15 les circuits croisillons se superposent a l'interieur 
de chaque bloc de calcul. 

Les circuits croisillons auxquels est affecte le 
coefficient 1/0, sont appeles croisillons de bord car 

20 ils effectuent des calculs sur les echantillons 
disposes aux extremites du bloc de calcul. Les autres 
circuits croisillons sont appeles circuits croisillons 
internes. II convient de noter, qu'a chaque etape de 
transformation, on ne traite pas tou jours la totalite 

25 des echantillons et que les echantillons non traites 
sont conserves dans leurs emplacements memoire pour 
etre traites a des etapes suivantes, ou bien produits 
en sortie s'ils sont deja dans leur forme finale. 

30 Le procede de transformation ainsi modifie delivrent 
des echantillons y(n) dans I'ordre croissant de 
l'indice n et presente a chaque etape de transformation 
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une symetrie de calcul facilitant l'adressage des 
echantillons a traiter. 

Selon un autre aspect de 1 1 invention, le procede de 
5 calcul comporte avantageusement m~1 etapes de 
transformation. Plusieurs modes de realisation derives 
de celui de la figure 8 et comprenant fi-l etapes de 
transformation sont done presentes dans la suite de la 
description, Un modele de circuit croisillon est 
10 associe a chacun de ces modes de realisation, 

Tous ces modes de realisation ont en commun les 
caracteristiques suivantes : 

- a chaque etape de transformation, il est prevu 
15 N/2 P+2 blocs de calcul et chaque bloc de calcul 

comporte un circuit croisillon de bord et/ou 2P-1 
circuits croisillons internes; tous les circuits 
croisillons, qu'ils soient de bord ou internes, 
effectuent des operations de calcul sur quatre 
20 echantillons reels; 

- les rangs des echantillons traites par un meme 
circuit croisillon sont definis de la maniere suivante: 
si on considere, a 1 ' etape de transformation E^ , un 
circuit croisillon de bord appartenant au bloc de 

25 calcul de rang a, il transforme les echantillons 
d'entree de rang 2^ +2 a, 2^ +2 a+2^ +1 -l , 2^ +2 a+2^ +1 , 
2^ +2 a+2^ +2 -l en des echantillons de sortie de meme rang 
et, si on considere un circuit croisillon interne de 
rang r dans le bloc de calcul de rang a a 1» etape E^, 

30 il transforme les echantillons d ' entree de rang 
2' ? - f2 a+2T+l, 2^ +2 a+2r+2, 2^ +2 a+2^ +2 -2r-3 , 2^ +2 a+2^ +2 -2T- 
2 en des echantillons de sortie de meme rang, avec /?>!; 
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- le coefficient affecte au circuit croisillon 
interne de rang r du bloc de calcul de rang a a I'etape 
Ep est egal a VJ 6 avec <S = (t+1) . (N/2^ +2 ) . 

5 Dans tous ces modes de realisation, les echantillons 
d' entree de chaque circuit croisillon forinent des 
couples d 1 echantillons, les rangs des echantillons d"un 
meme couple dans la suite des echantillons d 1 entree 
d'une etape de transformation etant symetriques par 

10 rapport a la valeur milieu des rangs extremes des 
echantillons d ' entree transformes par ledit circuit 
croisillon. Ce milieu correspond a la valeur 
2^ +2 a+2^ +1 -l/2. II suffit done de connaitre le rang de 
deux des quatres echantillons a appliquer sur les 

15 entrees du circuit croisillon pour en deduire le rang 
des deux autres. L'adressage de ces echantillons s'en 
trouve done simplifie. II sera detaille plus loin dans 
la description. 

20 Ainsi, selon un premier mode de realisation, les blocs 
de calcul voisins a chaque etape de transformation sont 
regroupes deux par deux, Les circuits croisillons de 
bord d'un meme bloc de calcul sont alors fusionnes en 
un circuit croisillon de bord unique. Un exemple de 

25 fusion de deux circuits croisillons de bord est 
presente a la figure 9. Cet exemple concerne les 
circuits croisillons de bord associes aux echantillons 
x(0) , x(8) , x(4) et x(12) . 

30 Par ailleurs, les etapes de transformation ne traitant 
pas a chaque fois la totalite des echantillons, 
certaines operations de calcul peuvent etre anticipees. 
Par exemple, le calcul des echantillons C R (l), C x (l) , 
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C R (3) et Cj(3) peut etre effectue a la deuxieme etape 
de transformation. II en resulte le schema de la figure 
10 presentant un premier mode de realisation du circuit 
de transformation dans laquelle le procede ne comporte 
5 que /z-1 etapes de transformation. A 1' etape de 
transformation E p/ chaque bloc de calcul comporte un 
circuit croisillon de bord et 2 p -l circuits croisillons 
internes. On peut considerer que ce procede comporte 
uniquement 3 etapes de transformation, la quatrieme 
10 etape se limitant a effectuer une addition et une 
soustraction. Cette addition et cette soustraction sont 
de preference effectuees lors de 1 • etape finale, et 
cela afin de limiter le nombre des etapes de 
transformation . 

15 

Un modele de circuit croisillon associe au mode de 
realisation de la figure 10 est represents a la figure 
11. II comporte: 

- quatre entrees pour recevoir des echantillons 
2 0 d' entree el, e2 , e3 , e4 et quatre sorties pour delivrer 

des echantillons de sortie si, s2, s3 , s4, et 

- trois entrees supplementaires respectivement de 
mode primaire MP, de permutation PERM et de coefficient 
COEF. 

25 

Ce circuit croisillon est charge d'appliquer 
selectivement sur les echantillons d 1 entree el, e2 , e3 
et e4, des operations de transformation differentes 
determinees chacune par les valeurs attributes a des 
30 signaux de mode primaire, de permutation et au 
coefficient W s admis sur les entrees supplementaires 
correspondantes . 
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Le signal de mode primaire vaut 0 s'il s'agit d'un 
circuit croisillon de bord et 1 s'il s'agit d'un 
circuit croisillon interne. Lorsque le signal de 
permutation est un 1, les echantillons de sortie si et 
5 s2 du circuit croisillon sont permutes avec les 
echantillons de sortie s3 et s4 . Cette permutation 
n'est possible que lorsque le circuit croisillon est un 
circuit croisillon interne, Enfin, le coefficient W s 
associe au circuit croisillon est applique sur 1' entree 
10 de coefficient COEF, 



Ainsi, si le coefficient complexe w s =A+j.B est applique 
sur 1 ' entree de coefficient du circuit croisillon, ce 
dernier delivre les echantillons de sortie si, s2 , s3 
15 et s4 suivants 

(1) si le signal de mode primaire est a 0: 

si = el + e2 
s2 = el - e2 
s3 = e4 - e3 

20 s4 = e3 + e4 

(2) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
S3 = el - A.e3 + B.e4 
25 s4 = -e2 + B.e3 + A.e4 

(3) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: sl=el-A.e3+B.e4 

s2 = -e2 + B.e3 + A.e4 
S3 = el + A.e3 -B.e4 
30 s4 = e2 + B.e3 + A.e4 



Selon un mode de realisation en variante, on peut 
prevoir d'effectuer l f addition et la soustraction dans 
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le circuit croisillon de bord de la derniere etape de 
transformation. Get mode de realisation est presente a 
la figure 12. A cet effet, le modele de circuit 
croisillon correspondant comporte une quatrieme entree 
5 supplementaire appelee entree de mode secondaire MS sur 
laquelle est applique un signal de mode secondaire. Ce 
signal est a 1 quand il s'agit du circuit croisillon de 
bord de la derniere etape de transformation, sinon il 
est a 0. Ce modele est illustre figure 13. Ce modele 

10 possede un mode de fonctionnement supplementaire par 
rapport au modele precedent; ainsi lorsque le signal de 
mode primaire est a 0 et que le signal de mode 
secondaire est a 1, on obtient en sortie: 
si = el + e2 + e3 + e4 

15 s2 = el - e2 

S3 = e4 - e3 

s4 = (el + e2) -(e3 + e4) 

Selon un deuxieme mode de realisation derive du schema 
20 de la figure 8, les circuits croisillons de bord 
voisins aux etapes de transformation E d'indice p pair 
sont regroupes deux par deux et sont fusionnes avec le 
circuit croisillon de bord de l 1 etape suivante d'indice 
impair pour former un nouveau circuit croisillon de 
25 bord au niveau de 1 1 etape d'indice impair. Ce 
regroupement est illustre a la figure 14 au travers 
d'un exemple. Dans cet exemple, les circuits 
croisillons de bord de la premiere etape de 
transformation traitant les echantillons x(0) , x(8), 
30 x(4) et x(12) sont fusionnes avec le circuit croisillon 
de bord de la deuxieme etape traitant les echantillons 
A R (0) et A R (2) . De meme, les circuits croisillons de 
bord de la premiere etape de transformation traitant 
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les echantillons x(2) , x(10) , x(6) et x(l4) sont 
fusionnes avec le circuit croisillon de bord de la 
deuxieme etape traitant les echantillons A R (4) et 
A R (6) . Les deux circuits croisi lions obtenus se 
dif f erencient en ce que le deuxieme effectue en plus 
une permutation entre la premiere et la derniere 
sortie. Si on applique ce regroupement a 1» ensemble de 
la transformation illustree figure 8, il en resulte que 
1' etape E Q n 1 utilise plus de circuits croisillons et 
peut etre supprimee. Le procede de transformation 
resultant est illustre a la figure 15. 



II convient cependant de distinguer deux cas pour ce 
procede de transformation : le cas ou N est une 
puissance paire de deux (/x pair) et le cas ou N est une 
puissance impaire de deux (n impair) . 

Dans le cas ou tx est pair, il y a un nombre pair 
d' etapes de transformation dans le mode de realisation 
presente a la figure 8 et le regroupement des circuits 
croisillons de bord des etapes d'indice pair avec ceux 
des etapes suivantes d'indice impair ne pose pas de 
probleme. Ce cas correspond au schema de la figure 15. 

Dans le cas oii /x est impair, le circuit croisillon de 
bord de la derniere etape d'indice pair ne peut etre 
regroupe avec d'autres circuits croisillons de bord. II 
convient alors de prevoir un mode de f onctionnement 
particulier pour ce cas. Ce cas est illustre figure 16, 
laquelle figure represente la transformation d'une 
suite de huit echantillons reels (fi=3) . Le circuit 
croisillon de bord de la derniere etape de 
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transformation de ce circuit n'a pu etre fusionne avec 
d 1 autres circuits croisillons de bord. 

Le modele de circuit croisillon associe a ce deuxieme 
5 mode de realisation est illustre figure 17; il se 
differencie du modele precedent en ce que le signal de 
mode secondaire est a 1 lorsqu'il s'agit d'un circuit 
croisillon de bord destine a mettre en oeuvre la 
derniere etape et que /x est impair, et en ce que la 
10 permutation s 1 applique a 1' ensemble des circuits 
croisillons des blocs de calcul de rang pair* 

Les operations de calcul realisees par le circuit 
croisillon sont egalement differentes et sont definies 
15 de la maniere suivante : 

(1) si les signaux de mode primaire, de mode secondaire 
et de permutation sont a 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 

20 s3 = e4 - e3 

s4 = (el + e2)-(e3 + e4) 

(2) si le signal de mode primaire est a 0 et le signal 
de mode secondaire est a 1: 

si = el + e4 
25 s2 = e2 

S3 = e3 
s4 = el - e4 

(3) si le signal de mode primaire est a 0 et le signal 
de permutation a l: si = (e3 + e4)-(el+e2) 

30 s2 = el - e2 

s3 = e4 - e3 
s4 = el + e2 + e3 + e4 



m 
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(4) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

(5) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 4- A.e4 
S3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A.e4 



Pour le cas ou /x est impair, on peut egalement prevoir 
d'une part, de regrouper les circuits croisillons de 
bord pour la mise en oeuvre de la premiere etape de 
15 transformation de la meme fa?on que dans le premier 
mode de realisation et d 1 autre part, de regrouper les 
circuits croisillons des autres etapes comme dans le 
troisieme mode de realisation. Le regroupement des 
circuits croisillons a partir de la seconde 6tape est 
alors effectue en attribuant un indice pair a la 
premiere etape de transformation. Ces regroupements 
sont representes a la figure 18. 



Le modele de circuit croisillon correspondant a ce mode 
de realisation est represente a la figure 19. Le signal 
de mode secondaire vaut l s'il s'agit d'un circuit 
croisillon de bord mettant en oeuvre la premiere etape 
de transformation du circuit et si js est pair. Les 
operations de calcul effectuees par ce circuit 
30 croisillon sont les suivantes : 

(1) si le signal de mode primaire est a 0 et le signal 
de mode secondaire est a l: 

si = el + e2 
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s2 = el - e2 
S3 = e4 - e3 
s4 = e3 + e4 

(2) si les signaux de mode primaire, de mode secondaire 
5 et de permutation sont a 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 
s4 = (el + e2)-(e3 + e4) 
10 (3) si les signaux de mode primaire et de mode 
secondaire sont a 0 et le signal de permutation a 1: 

si = (e3 + e4)-(el + e2) 
s2 = el - e2 
s3 = e4 - e3 

15 s4 = el + e2 + e3 + e4 

(4) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
S3 = el - A.e3 + B.e4 
20 s4 = -e2 + B.e3 + A.e4 

(5) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: si = el - A,e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
S3 = el + A.e3 -B.e4 
25 s4 = e2 + B.e3 + A.e4 



Dans tous les modes de realisation decrits 
precedemment, les rangs des echantillons traites par un 
meme circuit croisillon sont symetriques deux a deux 
30 par rapport a une valeur milieu. II suffit done de 
connaitre le rang des deux premiers echantillons 
d 1 entree du circuit croisillon pour en deduire les deux 
autres par symetrie. Si les echantillons d' entree et de 
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sortie sont sauvegardes dans des emplacements memoire 
dont 1 1 adresse correspond au rang de ces echantillons , 
l'adressage de ces derniers s'en trouve simplifie. 

5 En effet, il suffit alors de generer deux adresses par 
circuit croisillon, les deux autres pouvant etre 
deduites par symetrie. Par ailleurs, il convient de 
noter que les adresses des echantillons d' entree et 
celles des echantillons de sortie correspondants sont 

0 les memes car la transformation est appliquee selon une 
methode "in place". 



Les adresses associees aux differents circuits 
croisillons se rapportant au procede de transformation 

15 de la figure 12 sont presentes a la figure 20. 

L 1 adresse d'un echantillon est prise egale au rang de 
cet echantillon dans la suite des echantillons a 
laquelle il appartient. La suite a transformer dans 
l'exemple de la figure 12 comporte 16 echantillons. II 

20 est done necessaire de produire 16 adresses, 
l'adressage se faisant sur quatre bits. Pour une suite 
de N echantillons, l'adressage se fait sur log 2 (N) 
bits . 



Chaque etape de transformation du procede de la figure 
12 est mise en oeuvre par quatre circuits croisillons 
traitant chacun quatre echantillons reels d' entree. Les 
adresses binaires des echantillons a traiter a chaque 
etape sont done reparties en quatre groupes de quatre 
adresses. Les groupes d 1 adresses se rapportant a un 
circuit croisillon de bord sont contenus dans des 
rectangles en trait gras et les groupes d 1 adresses se 
rapportant a un circuit croisillon interne sont 
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contenus dans des rectangles en trait fin. Par 
ailleurs, les groupes d'adresses se rapportant a un 
meme bloc de calcul sont regroupes dans des rectangles 
en pointilles. 

5 

Le premier circuit croisillon de bord a l'etape E Q 
traite les guatre premiers echantillons de la suite 
x(n) issue de l'etape preliminaire de classement. Les 
adresses correspondantes a generer pour ce circuit 

10 croisillon sont done 0000, 0001, 0010 et 0011. De meme 
si on considere le premier circuit croisillon de bord 
de l'etape E x , il traite les echantillons de rang 0, 3, 
4 et 7 de la suite A(n). Les adresses correspondantes a 
generer pour ce circuit croisillon sont done 0000, 

15 0011, 0100 et 0111. 

Les (AX-2-p) bits de poids fort des adresses indiquent 
le rang q du bloc de calcul auquel se rapporte 
l'adresse. Ainsi, les deux bits de poids fort des 

2 0 adresses generees pour le premier bloc de calcul de 

l*etape E Q sont 00. II faut noter egalement que, etant 
donne que la derniere etape n' utilise qu'un seul bloc 
de calcul, les adresses generees pour ce bloc de calcul 
ne comporte pas de bit se rapportant au rang de ce bloc 
25 (/x-p-2=0) . 

Pour chaque circuit croisillon, on prevoir de generer 
uniquement deux adresses binaires; les deux autres sont 
obtenus en inversant les (p+2) bits de poids faible des 

3 0 adresses generees. Si on considere par exemple les 

adresses du premier circuit croisillon de bord a 
l'etape E lf on ne genere que les adresses 0000 et 0011 
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et on obtient les adresses 0111 et 0100 en inversant 
les 3 bits de poids faible des adresses generees . 

Une premiere et une seconde adresse par circuit 
croisillon sont alors produites, lesquelles adresses 
sont consecutives lorsqu'il s'agit d'un circuit 
croisillon interne. Lorsqu'il s'agit d'un circuit 
croisillon de bord, les p+2 bits de poids faible de la 
premiere adresse sont egaux a 0, et les p+2 bits de 
poids faible de la seconde adresse forment un nombre 
egal a 2P +1 -l. 



En ce qui concerne les coefficients W s a appliquer sur 
les entrees de coefficient COEF des circuits 
croisillons internes, ils sont stockes dans une memoire 
du circuit de calcul. Seules N/4 valeurs de 
coefficients sont necessaires au calcul de la 
transformee de Fourier. L'adressage de ces coefficients 
s'effectue sur bits. Dans les modes de realisation 

decrits precedemment ou N=16, on utilise uniquement les 
couples 1/0, 2/ -4, 3/ -3 et 4/ -2 correspondant aux 
coefficients W°, W 1 , W 2 , w 3 . 

Selon 1' invention, il est done necessaire de connaitre 
1« adresse de ces quatre coefficients. L' adresse 
associee a chaque coefficient W s est choisie egale a la 
valeur de la puissance s. Par consequent, les adresses 
des coefficients W°, W 1 , w 2 , W 3 sont respectivement 00, 
01, 10, 11. 

Afin de ne pas avoir a generer ces adresses, on utilise 
selon 1' invention les adresses produites pour adresser 
les echantillons. L 1 adresse du coefficient qui est 
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affecte a un circuit croisillon est comprise dans la 
seconde adresse produite correspondant a la plus elevee 
des deux adresses. 

5 Cependant, il convient de distinguer trois cas: 

(a) lorsque p+l=M-2, !• adresse du coefficient 
correspond au nombre forme par les p+1 bits de poids 
faible de la seconde adresse generee pour ce circuit 
croisillon interne. C'est le cas pour la deuxieme etape 

10 (P=l) dans I'exemple de la figure 20. Les 2 bits de 
poids faible de la seconde adresse sont 10 et designent 
done le coefficient VI 2 . 

(b) si p+l>/x-2, 1" adresse du coefficient correspond au 
nombre forme par les p+l bits de poids faible de la 

15 seconde adresse generee pour ce circuit croisillon 
interne, ote de ses ju-p-l bits de poids faible. C'est 
le cas pour la troisieme etape (p=2) dans l'exemple de 
la figure 20. Les 3 bits de poids faible de la seconde 
adresse generee pour le premier circuit croisillon 

2 0 interne sont 010. Lorsqu'on ote le dernier bit 
()Li-p-l=l) de ce nombre, on obtient le nombre 01 se 
rapportant au coefficient w 1 . Ce cas correspond 
tou jours a la dernier e etape d'un procede de 
transformation comportant jlz-1 etapes de transformation. 

25 (c) si p+ K/x-2 , 1' adresse du coefficient correspond au 
nombre forme par les p+1 bits de poids faible de la 
seconde adresse generee pour ce circuit croisillon 
interne, suivis de M-p-3 bits a zero en fin de nombre. 
Ce cas est illustre figure 21. Cette figure represente 

30 les adresses se rapportant a un circuit croisillon 
interne pour la mise en oeuvre de la deuxieme etape de 
transformation (p=l) d'un procede de transformation 
destine a traiter une suite de 3 2 echantillons reels 
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(N=32 et m=5) . Ce circuit croisillon transforme les 
echantillons de rang 1, 2, 5 et 6 de la suite des 
echantillons obtenue en fin de premiere etape. Les deux 
bits de poids faibles de la seconde adresse generee 
5 sont 10 et si on y rajoute un zero (M~p-3=1) en fin de 
nombre, on obtient le nombre 100 designant le 
coefficient W 4 * 

Ainsi, les deux adresses produites par le generateur 
10 d 1 adresses pour un circuit croisillon servent a 
adresser a la fois les quatre echantillons a traiter et 
le coefficient se rapportant au circuit croisillon. 

De preference, les echantillons d' adresse paire et 
15 d 1 adresse impaire seront stockes dans deux memoires 
distinctes. Ainsi, deux echantillons d 1 entree pourront 
etre lus simultanement et les echantillons de sortie 
resultants pourront etre ecrits simultanement, ce qui 
represente un gain au niveau du temps de traitement de 
20 la suite a transformer. 

Comme mentionne precedemment , a chaque procede de 
calcul fonctionnant selon un algorithme avec 
entrelacement temporel correspond un procede 

25 fonctionnant selon un algorithme avec entrelacement 
frequentiel. Pour l'obtenir, il suffit d'une part 
d'inverser la suite des operations de transformation du 
procede avec entrelacement temporel correspondant , et 
d 1 autre part, pour chaque circuit croisillon, 

30 d'inverser les operations de transformation en elles- 
memes par rapport a celle du procede correspondant. 
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Aussi, 1 1 invention concerne egalement un procede de 
calcul de la transformee de Fourier rapide ou de la 
transformee de Fourier rapide inverse d'une suite de N 
echantillons complexes X(n) conjugues deux a deux, avec 
5 N puissance de 2, fonctionnant selon un algorithme avec 
entrelacement frequentiel. La suite des N echantillons 
complexes X(n) est representee par une suite de N 
echantillons reels y(n) definie de la maniere suivante: 
y(0) = Re[X(0) ] 
10 y(n) = Re[X( (n+1) /2) ] pour n impair et 

different de N-l 
y(n) = Im[X(n/2)] pour n pair et different 

de 0 

y(N-l) = Re[X(N/2) ] 

15 

Selon 1' invention, ce procede comporte essentiellement 
des etapes de transformation pour transformer des 
echantillons d ' entree en echantillons de sortie. Les 
echantillons reels y(n) sont traites a une premiere 

2 0 etape de transformation et la derniere etape delivre 
une suite de N echantillons de sortie x(n) 
representatif s de la transformee de Fourier rapide ou 
rapide inverse de la suite des echantillons X(n) . Comme 
pour les circuits avec entrelacement temporel, chaque 

25 etape de transformation est raise en oeuvre par un jeu 
de circuits croisillons a plusieurs entrees et a 
plusieurs sorties. La transformation etant effectuee 
selon une methode "in place 11 , 1' ensemble des etapes est 
effectue au moyen d'un unique jeu de circuits 

30 croisillons dont le mode operatoire est modifie a 
chaque etape de transformation. Les echantillons 
d 1 entree et de sortie sont memorises a chaque etape de 
transformation dans une memoire de stockage. Apres 
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transformation, les echantillons de sortie d'un meme 
circuit croisillon se substituent dans la memoire de 
stockage aux echantillons d' entree de meme rang 
correspondants . 

5 

Selon 1' invention, si les echantillons y(n) presentes a 
la premiere etape de transformation sont classes dans 
l'ordre croissant de 1 • indice n, les echantillons de 
sortie x(n) sont delivres a la derniere etape de 
10 transformation dans l'ordre binaire inverse de 1 1 indice 
n. L'ordre de sortie des echantillons x(n) peut ensuite 
etre modifie par une etape finale de maniere a les 
classer dans l'ordre croissant de 1' indice n. 

Un mode de realisation d'une telle transformation est 
represents a la figure 22. II se deduit du mode de 
realisation de la figure 12 en inversant I'agencement 
fonctionnel des etapes de transformation de la figure 
12 (image en miroir de ce qu'il est pour un 
entrelacement temporel) • Ce mode de realisation permet 
de calculer la transformee de Fourier rapide inverse 
d'une suite de 16 echantillons reels y(n) 
representatif s d'une suite de 16 echantillons complexes 
X(n) conjugues deux a deux. Les coefficients W s sont 
par consequent du type e^ ( 27rs /N) # 

Ce mode de realisation comprend trois etapes de 
transformation E p avec 0<p<2 . A chaque etape de 
transformation E p/ les circuits croisillons sont 
30 desormais repartis en 2? blocs de calcul, ces blocs de 
calcul etant ordonnes a chaque etape selon un rang q 
croissant allant de 0 a 2 P -1. Chaque bloc de calcul 
comporte un circuit croisillon de bord et N/2 p+2 -l 
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circuits croisillons internes, Les echantillons y(n) 
sont appliques dans l*ordre de l^ndice n a la premiere 
etape. 

5 Les rangs des echantillons traites par un meme circuit 
croisillon sont definis de la maniere suivante: si on 
considere a 1* etape un circuit croisillon de bord 

appartenant au bloc de calcul de rang a, il transforme 
les echantillons d 1 entree de rang 2^"^ a, 2* x ~^a+2* i ~^~ 1 - 

10 1, 2^"^ 3 a+2^-/ 3 - 1 / 2^a+2^-l en des echantillons de 
sortie de meme rang et, si on considere a 1 1 etape E^ un 
circuit croisillon interne de rang r dans le bloc de 
calcul de rang a, il transforme les echantillons 
d 1 entree de rang 2* i ~0a+2r+i / 2 M ~^a+2r+2 / 2^~^a+2^"^-2T- 

15 3, 2^~^a+2 /x ~^-2r-2 en des echantillons de sortie de 
meme rang. Enfin, le coefficient affecte au circuit 
croisillon interne de rang t du bloc de calcul de rang 
a a l'etape E^ est egal a VI s avec 6 = (t+1).2^. 

2 0 Les echantillons reels x(n) obtenus en fin de procede 
sont delivres dans 1 ' ordre binaire inverse de l'indice 
n . 

Au niveau des circuits croisillons, les coefficients W s 
25 sont du type e 3 < 27rs / N ) et, les entrees et les sorties 
ont ete interverties par rapport au mode de realisation 
de la figure 12. En consequence, les operations 
effectuees par les circuits croisillons de ce mode de 
realisation sont differentes de celles effectuees par 
30 les circuits croisillons de la figure 12. Un modele de 
circuit croisillon associe au mode de realisation de la 
figure 22 est represents a la figure 23. 
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Tout comme le modele associe a la figure 12, il 
comporte quatre entrees et quatre sorties de donnees 
ainsi que quatre entrees supplementaires respectivement 
de mode primaire MP, de mode secondaire MS, de 
5 permutation PERM et de coefficient COEF. Le signal de 
mode primaire vaut 0 s'il s'agit d'un circuit 
croisillon de bord et 1 s'il s'agit d'un circuit 
croisillon interne, Le signal de permutation vaut 0 
pour les valeurs paires du rang q et 1 pour les valeurs 
10 impaires. Enfin, le signal de mode secondaire vaut 1 si 
le circuit croisillon de bord est utilise pour la mise 
en oeuvre de la premiere etape et 0 sinon. 

Les operations de calcul effectuees par ce circuit 
15 croisillon sont les suivantes (W s =A+j.B): 

(1) si les signaux de mode primaire et de mode 
secondaire sont a 0: si = (el + e2)/2 

s2 = (el - e2) /2 
S3 = (e4 - e3) /2 
20 s4 = (e3 + e4) /2 

(2) si le signal de mode primaire est a 0 et le signal 
de mode secondaire est a 1: si = [ (el+ e4 ) /2+e2 ] /2 

s2 = [ (el+e4) /2-e2] /2 
S3 = -[e3-(el-e4)/2]/2 
25 s4 = [e3+(el-e4) /2] /2 

(3) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si =(el+e3) /2 

s2 =(e2-e4)/2 

S3 =[(el-e3) .A - (e2+e4).B]/2 
30 s4 =[-(el-e3).B + (e2+e4 ) . A] /2 

(4) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: si =[ (el-e3) .A - (e2+e4).B]/2 

s2 =[-(el-e3).B +(e2-fe4) . A] /2 
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S3 =(el+e3)/2 
s4 =(e2-e4)/2 

Ce modele est deduit du modele de la figure 13 en 
5 intervertissant les valeurs el, e2 , e3, e4 et les 
valeurs si, s2 , s3, s4, et en remplagant B par -B car 
le coefficient W s est a present du type e3( 27rs / N ). 
Etant donne que A et B representent le cosinus et le 
sinus d'un itieme nombre, on a en fait A 2 +B 2 =l. Les 
10 expressions de si, s2, s3 , s4 s'en trouvent 
simplif iees . 

Selon I 1 invention, des adresses sont par ailleurs 
generees pour l'adressage des echantillons. II est 

15 prevu de generer deux adresses binaires de fi bits par 
circuit croisillon, chaque adresse binaire 

correspondant au rang d'un echantillon d ' entree du 
circuit croisillon. Les adresses des deux autres 
echantillons a appliquer sur les entrees du circuit 

20 croisillon sont obtenus en inversant les ju-p bits de 
poids faible des deux premieres adresses. 

De la meme fagon que pour les procedes de 
transformation fonctionnant selon un entrelacement 

25 temporel, les deux adresses binaires produites sont 
consecutives s*il s'agit d'un circuit croisillon 
interne. S'il s'agit d'un circuit croisillon de bord, 
les ju-p bits de poids faible de la premiere adresse 
generee sont egaux a 0, et les jiz-p bits de poids faible 

3 0 de la seconde adresse forment un nombre egal a 
N/2 P* hl -l. A titre d 1 illustration, les adresses 
produites pour le circuit de transformation de la 
figure 22 sont rassemblees a la figure 24. 
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Avantageusement, on pourra prevoir de stocker les 
echantillons d'adresse paire et les echantillons 
d'adresse impaire dans deux memoires distinctes afin de 
5 diminuer le temps de traitement de 1» operation de 
transformation . 

Enfin, les adresses generees pour adresser les 
echantillons sont egalement utilises pour adresser les 

10 coefficients W s . La valeur du parametre s est utilise 
pour adresser le coefficient W s correspondant . Dans cet 
mode de realisation, le parametre s est egal: 

- si M-p-l^M-2, au nombre forme par les /x-p-1 bits 
de poids faible de la seconde adresse produite pour 

15 ledit circuit croisillon interne, 

~ si /z-p-KjLt-2, au nombre forme par les ju-p-1 bits 
de poids faible de la seconde adresse produite pour 
ledit circuit croisillon interne, suivi de p-l bits a 
zero en fin de nombre, 

20 " si /i-p-l>/x-2, au nombre forme par les ju-p-1 bits 

de poids faible de la seconde adresse produite pour 
ledit circuit croisillon interne, ote de ses p+1 bits 
de poids faible. Ce cas correspond a la premiere etape 
(p=0) des procedes de transformation fonctionnant selon 

25 un algorithme avec entrelacement frequentiel. 
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R E VEND I CAT I ONS 

1. Procede de calcul de la transformee de Fourier 
rapide ou de la transformee de Fourier rapide inverse 
5 d f un signal numerique defini par une suite de N 
echantillons reels de depart x(n) , avec N puissance de 
deux et n e [0.*N-1], comportant des etapes successives 
de transformation (2) pour transformer des echantillons 
d' entree en echantillons de sortie, 1" ensemble des 

10 etapes de transformation etant effectue au moyen d'un 
unique jeu de circuits croisillons a plusieurs entrees 
et plusieurs sorties dont le mode operatoire est change 
select ivement a chaque etape de transformation, les 
echantillons d 1 entree et de sortie de chaque Stape de 

15 transformation etant stockes dans une memoire de 
stockage, une suite de N echantillons de sortie y (n) 
representatif s de la transformee de Fourier rapide ou 
rapide inverse des echantillons de depart x(n) etant 
delivree a la derniere etape de transformation, 

20 caracterise en ce que les echantillons de sortie 

y(n) sont reels, 

et en ce que les echantillons de sortie d'un 
circuit croisillon se substituent dans la memoire de 
stockage aux echantillons d' entree de meme rang 

25 correspondants, de telle sorte que, si les echantillons 
de depart x(n) traites a la premiere etape de 
transformation sont classes dans l f ordre binaire 
inverse de leur indice n, les echantillons de sortie 
y(n) sont delivres a la derniere etape de 

30 transformation dans 1 • ordre croissant de I 1 indice n, 
lesquels echantillons de sortie sont definis par les 
relations suivantes : 
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y(0) = Re[X(0)] 

y(n) = Re[X((n+l)/2) ] pour n impair et 

different de N-l 
y(n) = Im[X(n/2)] pour n pair et different 

5 de 0 

y(N-l) = Re[X(N/2) ] 
ou les echantillons X(n) , avec n e [0..N-1], designent 
les echantillons complexes de la suite correspondant a 
la transformee de Fourier rapide ou rapide inverse de 
10 la suite des echantillons de depart x(n). 

2. Procede de calcul de la transformee de Fourier 
rapide ou de la transformee de Fourier rapide inverse 
d'un signal numerique defini par une suite de N 
15 echantillons complexes X(n) conjugues deux a deux 
representes par une suite de N echantillons reels de 
depart y(n), avec N puissance de deux et n e [0..N-1], 
les Echantillons de depart y(n) etant definis de la 
maniere suivante : 
20 y(0) = Re[X(0) ] 

y(n) = Re[X((n+l)/2) ] pour n impair et 

different de N-l 
y(n) = Im[X(n/2)] pour n pair et different 

de 0 

25 y(N-l) = Re[X(N/2)] 

lequel procede de calcul comporte des etapes 
successives de transformation pour transformer des 
echantillons d« entree en echantillons de sortie, une 
suite de N echantillons de sortie x(n) representatif s 

3 0 de cette transformee de Fourier rapide ou rapide 
inverse etant delivree a la derniere etape de 
transformation, 1« ensemble des etapes de transformation 
etant effectue au moyen d'un unique jeu de circuits 
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croisillons a plusieurs entrees et plusieurs sorties 
dont le mode operatoire est change selectivement a 
chague etape de transformation, les echantillons 
d' entree et de sortie de chague etape de transformation 
5 etant stockes dans une memo ire de stockage, 

caracterise en ce gue les echantillons de sortie 
x(n) sont reels, 

et en ce gue les echantillons de sortie d'un 
circuit croisillon se substituent dans la memoire de 

10 stockage aux echantillons d 1 entree de meme rang 
correspondants, de telle sorte gue, si les echantillons 
de depart y(n) traites a la premiere etape de 
transformation sont classes dans 1 1 ordre croissant de 
1 ' indice n, les echantillons de sortie x(n) sont 

15 delivres a la derniere etape de transformation dans 
1 1 ordre binaire inverse de l 1 indice n. 

3. Procede de calcul selon la revendication 1 ou 2 , 
caracterise en ce gue, a chague etape de 

20 transformation; chague circuit croisillon transforme 
des couples d f echantillons d" entree, les rangs des 
echantillons d 1 entree d'un meme couple dans la suite 
des echantillons d* entree de ladite etape de 
transformation etant symetrigues par rapport a un 

25 milieu entre les valeurs de rang extremes des 
echantillons d 1 entree transformes par ledit circuit 
croisillon . 

4. Procede de calcul selon la revendication 3, 
30 caracterise en ce gu'il comporte /z-l etapes de 

transformation E avec /x=log 2 (N) et p € [0../X-2]. 
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5. Procede de calcul selon la revendication 4 elle-meme 
dependante de la revendication 3 elle-meme dependante 
de la revendication 1, caracterise en ce qu'il comporte 
en outre : 

5 - une etape preliminaire pour modifier l^rdre des 

echantillons de depart x(n) ranges dans l f ordre 
croissant de 1 1 indice n et les presenter dans l'ordre 
binaire inverse de 1» indice n a la premiere etape de 
trans f ormat ion , et 
10 - u ne etape finale pour traiter la suite 

d f echantillons de sortie y(n) et delivrer une suite de 
N echantillons complexes X(n) conjugues correspondant a 
la transformee de Fourier rapide ou rapide inverse de 
la suite des echantillons de depart x(n) . 

15 

6. Procede de calcul selon la revendication 4 elle-meme 
dependante de la revendication 3 elle-meme dependante 
de la revendication 1, ou selon la revendication 5, 
caracterise en ce que, a chaque etape de transformation 

2 0 E p , les circuits croisillons sont repartis en N/2P +2 
blocs de calcul, 

en ce que chaque bloc de calcul comporte un 
circuit croisillon de bord et/ou 2P-1 circuits 
croisillons internes , 

25 en ce que le circuit croisillon de bord du bloc de 

calcul de rang a a 1 1 etape de transformation 
transforme les echantillons d' entree de rang 2^ +2 a, 
S^a-K^-l, 2* + 2 a +2* + l, 20 +2 a + 20 +2 -l en des 
echantillons de sortie de meme rang, 

30 et en ce que le circuit croisillon interne de rang 

t du bloc de calcul de rang a a 1 ■ etape de 
transformation E^ transforme les echantillons d' entree 
de rang 2^ +2 a+2r+l , 2^ +2 a+2r+2 , 2^ +2 a+2^ +2 -2r-3 , 
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2 p " r ^a+2 p " r-£ -2r-2 en des echantillons de sortie de meme 
rang, avec /3>1. 

7. Precede de calcul selon la revendication 6, 
5 caracterise en ce que, a chaque circuit croisillon est 

affecte un coefficient W s , coefficient sur lequel est 
base 1' operation de calcul a 1' inter ieur du circuit 
croisillon, ledit coefficient etant egal a e~3 (27rs/N) 
avec s € [0.. N/4-1] dans le cas d'une transformee de 
10 Fourier rapide et est egal a e 3( 27rs / N ) avec s e 
[0.. N/4-1] dans le cas d'une transformee de Fourier 
rapide inverse. 

8. Procede de calcul selon la revendication 7, 
15 caracterise en ce que, au circuit croisillon interne de 

rang t du bloc de calcul de rang a a l'etape de 
transformation Eg est affecte le coefficient avec 
S = (t+1) . (N/2^ +2 ) . 

20 9, Procede de calcul selon la revendication 8, 
caracterise en ce que les circuits croisillons destines 
a mettre en oeuvre les etapes de transformation sont 
tous de meme type et comportent 

- quatre entrees pour recevoir des echantillons 
25 d 1 entree et quatre sorties pour delivrer des 

echantillons de sortie, 

- quatre entrees supplementaires respect ivement de 
mode primaire, de mode secondaire, de permutation et de 
coefficient, 

3 0 afin d'appliquer selectivement sur les 

echantillons d 1 entree, des operations de transformation 
differentes determinees chacune par les valeurs 
attributes a des signaux de mode primaire, de mode 
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secondaire, de permutation et a un coefficient admis 
sur lesdites entrees supplementaires correspondantes . 

10. Procede de calcul selon la revendication 9, 
caracterise en ce que, pour chaque circuit croisillon, 
le signal de mode primaire vaut 0 s'il s'agit d'un 
circuit croisillon de bord et 1 s'il s'agit d'un 
circuit croisillon interne, 

en ce que le signal de permutation vaut 0 pour les 
blocs de calcul de rang pair, y compris le rang 0, et 1 
pour les autres . 



11. Procede de calcul selon la revendication 10, 
caracterise en ce que, a I'etape de transformation Ep, 
chaque bloc de calcul comporte un circuit croisillon de 
bord et 2?-l circuits croisillons internes. 



12. Procede de calcul selon la revendication 11, 
caracterise en ce que le signal de mode secondaire vaut 
1 si le circuit croisillon de bord est utilise pour la 
dernier e etape de transformation et 0 sinon. 

13. Procede de calcul selon la revendication 12, 
caracterise en ce que, pour quatre echantillons 
d' entree el, e2 , e3 et e4 et pour un coefficient 
complexe W s =A+j.B, le circuit croisillon fournit les 
echantillons de sortie si, s2, s3 et s4 suivants 

(1) si les signaux de mode primaire et de mode 
secondaire sont a 0: si = el + e2 

s2 = el - e2 

S3 = e4 - e3 

s4 = e3 + e4 



• 
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(2) si le signal de mode primaire est a 0 et le signal 
de mode secondaire est a 1: si = el + e2 + e3 + e4 

s2 = el - e2 
s3 = e4 - e3 

5 s4 = (el + "e2) -(e3 + e4) 

(3) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
S3 = el - A.e3 + B.e4 
10 s4 = -e2 + B.e3 + A.e4 

(4) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 -B.e4 
15 s4 = e2 + B.e3 + A.e4 

14* Precede de calcul selon la revendication 10, 
caracterise en ce que, a 1 1 etape de transformation E , 
chaque bloc de calcul comporte: 
20 - 2 p -l circuits croisillons internes et un circuit 

croisillon de bord pour les valeurs paires de 1 1 indice 
p ainsi que pour la derniere etape de transformation si 
fj. est pair, et 

- 2 p -l circuits croisillons internes sinon. 

25 

15. Procede de calcul selon la revendication 13, 
caracterise en ce que, le signal de mode secondaire 
vaut 1 si le circuit croisillon de bord est utilise 
pour la derniere etape de transformation avec fi impair 

30 et 0 sinon. 

16. Procede de calcul selon la revendication 15, 
caracterise en ce que, pour quatre echantillons 
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d' entree el, e2, e3 et e4 et pour un coefficient 
complexe w s =A+j.B, le circuit croisillon fournit les 
echantillons de sortie si, s2, s3 et s4 suivants 

(1) si les signaux de mode primaire, de mode secondaire 
5 et de permutation sont a 0: si = el + e2 + e3 + e4 

s2 = el - e2 
s3 = e4 - e3 
s4 = (el + e2)-(e3 + e4) 

(2) si le signal de mode primaire est a 0 et le signal 
10 de mode secondaire est a 1: si = el + e4 

s2 = e2 
s3 = e3 
s4 = el - e4 

(3) si le signal de mode primaire est a 0 et le signal 
15 de permutation a 1: si = (e3 + e4)-(el+e2) 

s2 = el - e2 
s3 = e4 - e3 
s4 = el + e2 + e3 + e4 

(4) si le signal de mode primaire est a 1 et le signal 
20 de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
s3 = el - A,e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

(5) si le signal de mode primaire est a 1 et le signal 
2 5 de permutation a 1: si = el - A.e3 + B.e4 

S2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A.e4 



30 



17. Procede de calcul selon la revendication 10, 
caracterise en ce que, a l'etape de transformation E p , 
chaque bloc de calcul comporte: 
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- 2 p -l circuits croisillons internes et un circuit 
croisillon de bord pour les valeurs paires de l'indice 
P, et 

- 2 p -l circuits croisillons internes sinon. 

5 

18. Procede de calcul selon la revendication 17 , 
caracterise en ce que le signal de mode secondaire vaut 
1 si le circuit croisillon de bord est utilise pour la 
premiere etape de transformation avec /jl pair et 0 

10 sinon. 

19. Procede de calcul selon la revendication 18, 
caracterise en ce que, pour quatre echanti lions 
d 1 entree el, e2 , e3 et e4 et pour un coefficient 

15 complexe W s =A+j.B, le circuit croisillon fournit les 
echantillons de sortie si, s2 , s3 et s4 suivants 

(1) si le signal de mode primaire est a 0 et le signal 
de mode secondaire est a 1: si = el + e2 

s2 = el - e2 

20 s3 = e4 - e3 

s4 = e3 + e4 

(2) si les signaux de mode primaire, de mode secondaire 
et de permutation sont a 0: si = el + e2 + e3 + e4 

s2 = el - e2 

25 s3 = e4 - e3 

s4 = (el + e2)-(e3 + e4) 

(3) si les signaux de mode primaire et de mode 
secondaire sont a 0 et le signal de permutation a 1: 

si = (e3 + e4)-(el + e2) 
30 s2 = el - e2 

s3 = e4 - e3 
s4 = el + e2 + e3 + e4 
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(4) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
S3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

(5) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
S3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A.e4 



20. Procede de calcul selon la revendication 8, 
caracterise en ce que les circuits croisillons destines 
a mettre en oeuvre les etapes de transformation sont 
15 tous de meme type et comportent 

- quatre entrees pour recevoir des echantillons 
d« entree et quatre sorties pour delivrer des 
echantillons de sortie, 

- trois entrees supplementaires respectivement de 
2 0 mode primaire, de permutation et de coefficient, 

afin d'appliquer selectivement sur les 
echantillons d' entree, des operations de transformation 
differentes determinees chacune par les valeurs 
attributes a des signaux de mode primaire, de 

2 5 permutation et a un coefficient admis sur lesdites 

entrees supplementaires correspondantes , 

et en ce que l'etape finale effectue en outre une 
addition et une soustraction entre le premier et le 
dernier echantillon de sortie delivres a la derniere 

3 0 etape de transformation. 



21. Procede de calcul selon la revendication 20, 
caracterise en ce que, a !• etape de transformation E , 
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chaque bloc de calcul comporte un circuit croisillon de 
bord et 2 p -l circuits croisillons internes. 



22. Procede de calcul selon la revendication 21, 
5 caracterise en ce que, pour quatre echantillons 
d' entree el, e2 , e3 et e4 et pour un coefficient 
complexe W s =A+j.B, le circuit croisillon fournit les 
echantillons de sortie si, s2 , s3 et s4 suivants 

(1) si le signal de mode primaire est a 0: 
10 si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 

(2) si le signal de mode primaire est a 1 et le signal 
15 de permutation a 0: si = el + A.e3 -B.e4 

s2 = e2 + B.e3 + A.e4 
S3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

(3) si le signal de mode primaire est a 1 et le signal 
2 0 de permutation a 1: si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A.e4 



25 23. Procede de calcul selon la revendication 9 ou 20, 
caracterise en ce que des premiere et seconde adresses 
binaires de y, bits sont generees pour chaque circuit 
croisillon, chaque adresse binaire correspondant au 
rang d'un echantillon d' entree dudit circuit croisillon 

30 et la seconde adresse binaire etant superieure a la 
premiere adresse binaire. 
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24. Procede de calcul selon la revendication 23, 
caracterise en ce que lesdites premiere et seconde 
adresses binaires sont consecutives s J il s'agit d'un 
circuit croisillon interne. 

25. Procede de calcul selon la revendication 23 ou 24, 
caracterise en ce que, s'il s f agit d'un circuit 
croisillon de bord, les p+2 bits de poids faible de la 
premiere adresse sont egaux a 0, et les p+2 bits de 
poids faible de la seconde adresse forment un nombre 
egal a 2P +1 -1. 



26. Procede de calcul selon la revendication 24 ou 25, 
caracterise en ce que 1« adresse des deux autres 
echantillons a appliquer sur les entrees du circuit 
croisillon, qu'il soit de bord ou interne, sont obtenus 
en inversant les (p+2) bits de poids faible desdites 
premiere et seconde adresses produites. 

27. Procede de calcul selon la revendication 26, 
caracterise en ce que les echantillons d' adresse paire 
et les echantillons d' adresse impaire sont stockes dans 
deux memoires distinctes. 



28. Procede de calcul selon la revendication 25, 
caracterise en ce que la valeur du parametre s du 
coefficient W s affecte a un circuit croisillon interne 
a l'etape de transformation E p est code sur /u-2 bits, 
et correspond: 

- si p+l=/x-2, au nombre forme par les p+l bits de 
poids faible de la seconde adresse binaire produite 
pour ledit circuit croisillon interne, 
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- si p+K/i-2, au nombre forme par les p+l bits de 
poids faible de la seconde adresse binaire produite 
pour ledit circuit croisillon interne, suivi de M~p-3 
bits a zero en fin de nombre, 
5 - si p+l>/i-2, au nombre forme par les p+l bits de 

poids faible de la seconde adresse binaire produite 
pour ledit circuit croisillon interne, ote de ses M"P-1 
bits de poids faible. 



10 29 . Procede de calcul selon la revendication 4 elle- 
meme dependante de la revendication 3 elle-meme 
dependante de la revendication 2, caracterise en ce 
que, a chaque etape de transformation E p , les circuits 
croisillons sont repartis en 2 p blocs de calcul, 

15 en ce que chaque bloc de calcul comporte un 

circuit croisillon de bord et N/2 p+2 -l circuits 
croisillons internes , 

en ce que le circuit croisillon de bord du bloc de 
calcul de rang a k 1« etape de transformation 

20 transforme les echantillons d 1 entree de rang 2^"^a, 
2^a+2^- 1 -l , 2^-^a+2^" 1 , 2^a+2^-l en des 
echantillons de sortie de meme rang, 

et en ce que le circuit croisillon interne de rang 
r du bloc de calcul de rang a a l 1 etape de 

25 transformation E^ transforme les echantillons d' entree 
de rang 2^~^G£+2t+1, 2^"^a+2r+2, 2^"^a+2^~^-2r-3 , 
2^"^a+2^"^-2r-2 en des echantillons de sortie de meme 
rang. 



30 



30. Procede de calcul selon la revendication 29, 
caracterise en ce qu'il comporte en outre une etape 
finale pour modifier 1 • ordre des echantillons de sortie 
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delivres a la derniere etape de transformation et les 
classer dans l'ordre croissant de 1 » indice n. 

31. Procede de calcul selon la revendication 29 ou 30, 
5 caracterise en ce que, a chaque circuit croisillon est 

affecte un coefficient W s , coefficient sur lequel est 
base 1" operation de calcul a I'interieur du circuit 
croisillon, ledit coefficient etant egal a e"3 (27rs/N) 
avec s e [0.. N/4-1] dans le cas d'une transformee de 
10 Fourier rapide et est egal a € >J(2ns/U) avec s e 
[0..N/4-1] dans le cas d'une transformee de Fourier 
rapide inverse . 

32. Procede de calcul selon la revendication 31, 
15 caracterise en ce que, au circuit croisillon interne de 

rang r du bloc de calcul de rang a a I 1 etape de 
transformation est associe le coefficient avec 

S = (r+1) .2^. 



33. Procede de calcul selon la revendication 32, 
caracterise en ce que les circuits croisillons pour 
mettre en oeuvre les etapes de transformation sont tous 
de meme type et comportent 

- quatre entrees pour recevoir des echantillons 
d 1 entree et quatre sorties pour delivrer des 
echantillons de sortie, 

- quatre entrees supplementaires respectivement de 
mode primaire, de mode secondaire, de permutation et de 
coefficient, 

afin d'appliquer selectivement sur les 
echantillons d f entree, des operations de transformation 
differentes determinees chacune par les valeurs 
attributes a des signaux de mode primaire, de mode 
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secondaire, de permutation et a un coefficient admis 
sur lesdites entrees supplementaires correspondantes . 

34. Procede de calcul selon la revendication 33, 
5 caracterise en ce que, pour chaque circuit croisillon, 

le signal de mode primaire vaut 0 s'il s'agit d'un 
circuit croisillon de bord et 1 s'il s'agit d'un 
circuit croisillon interne, 

et en ce que le signal de permutation vaut 0 pour 
10 les blocs de calcul de rang pair, y compris le rang 0, 
et 1 pour les valeurs impaires. 

35. Procede de calcul selon la revendication 31 ou 34, 
caracterise en ce que le signal de mode secondaire vaut 

15 1 si le circuit croisillon, qu'il soit de bord ou 
interne, est utilise pour la premiere etape de 
transformation et 0 sinon. 

36. Procede de calcul selon la revendication 35, 
2 0 caracterise en ce que, pour quatre echantillons 

d' entree el, e2, e3 et e4 et pour un coefficient 
complexe W s =A+j.B, le circuit croisillon fournit les 
echantillons de sortie si, s2, s3 et s4 suivants 

(1) si les signaux de mode primaire et de mode 
25 secondaire sont a 0: si = (el + e2) /2 

s2 = (el - e2)/2 
S3 = (e4 - e3)/2 
s4 = (e3 + e4) /2 

(2) si le signal de mode primaire est a 0 et le signal 
30 de mode secondaire est a 1: si = [ (el+e4 ) /2+e2 ] /2 

s2 = [ (el+e4) /2-e2]/2 
S3 = ~[e3-(el-e4) /2]/2 
s4 = [e3+(el-e4)/2]/2 
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(3) si le signal de mode primaire est a 1 et le signal 
de permutation a 0: si =(el-f-e3) /2 

s2 =(e2-e4)/2 

s3 =[(el-e3).A - (e2+e4).B]/2 
5 s4 =[-(el-e3).B +(e2+e4) .A]/2 

(4) si le signal de mode primaire est a 1 et le signal 
de permutation a 1: si =[(el-e3).A - (e2+e4).B]/2 

s2 =[-(el-e3).B +(e2+e4) .A]/2 
S3 =(el+e3)/2 

10 s4 =(e2-e4) /2 

37. Procede de calcul selon la revendication 33 , 
caracterise en ce que des premiere et seconde adresses 
binaires de /x bits sont generees pour chague circuit 
15 croisillon, chaque adresse binaire correspondant au 
rang d'un echantillon d 1 entree dudit circuit croisillon 
et la seconde adresse binaire etant superieure a la 
premiere adresse binaire. 

20 38. Procede de calcul selon la revendication 37, 
caracterise en ce que lesdites premiere et seconde 
adresses binaires sont consecutives s'il s'agit d'un 
circuit croisillon interne. 

25 39. Procede de calcul selon la revendication 37 ou 38, 
caracterise en ce que, s'il s'agit d'un circuit 
croisillon de bord, les /i-p bits de poids faible de la 
premiere adresse sont egaux a 0, et les /z-p bits de 
poids faible de la seconde adresse forment un nombre 

30 egal a N/2P +1 -1. 



40. Procede de calcul selon la revendication 38 ou 39, 
caracterise en ce que 1 ' adresse des deux autres 
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echantillons a appliguer sur les entrees du circuit 
croisillon sont obtenus en inversant les /i-p bits de 
poids faible des deux adresses produites. 

5 41. Precede de calcul selon la revendication 40, 
caracterise en ce que les echantillons d' adresse paire 
et les echantillons d' adresse impaire sont stockes dans 
deux memoires distinctes. 

10 42. Procede de calcul selon la revendication 41, 
caracterise en ce que la valeur du parametre s du 
coefficient W s affecte a un circuit croisillon interne 
a l'etape de transformation E est code sur bits, 
et correspond: 

15 - si Ai-p-l-/i-2, au nombre forme par les /x-p-l bits 

de poids faible de la seconde adresse produite pour 
ledit circuit croisillon interne, 

- si au nombre forme par les bits 
de poids faible de la seconde adresse produite pour 

20 ledit circuit croisillon interne, suivi de p-1 bits a 
zero en fin de nombre, 

- si jx-p-l>/x-2, au nombre forme par les |i-p-l bits 
de poids faible de la seconde adresse produite pour 
ledit circuit croisillon interne, ote de ses p+1 bits 

25 de poids faible. 
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[x(0) ,X(1) ,x(2) ,x(3) ,x(4) ,x(5) ,X(6) ,x(7) , 
X(8) ,X(9) ,X(10) ,x(ll) ,x(12) ,x(13) ,x(14) ,x(15)] 



> 




ETAPE PRELIMINAIRE I 


> 


f 



[x(0) ,x(8) ,x(4) ,x(12) ,x(2) ,x(10) ,x(6) ,x(14) , 
X(l) ,x(9) ,x(5) ,X(13) ,x(3) ,X(11) ,x(7) ,x(15) ] 







ETAPE S DE TRANSFORMATION 


> 





[y(O) ,y(l) ,y(2) ,y(3) ,y(4) ,y(5) ,y(6) ,y(7) , 
y(8) ,y(9) ,y(iO) ,y(ll) ,y(i2) ,y(i3) /y (i4) ,y(i5) ] 



> 




ETAPE 


FINALE 



[X(O) ,X(1) ,X(2) ,X(3) ,X(4) ,X(5) ,X(6) ,X(7) , 
X(8),X(9),X(10),X(11),X(12) f X(13) ,X(14) ,X(15)] 
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A R (4)+Aw(6)=B«(4) 
A w (4)-A w (6)=Br(6) 




•tA a (4)-A w (6))=-B«(6) 



Ap(4)+A,(6)=B«(4) 



FIG.7A 




r B R (9)^3.B R (13)>£B7(13)=C R (ir ~ 



B,(9).3.8 R (13)+3.B I (13)=C I Q) 



B R (9)-3.B R (13)-3.B,(13)=C R (11) 
-B,(9)-3,B R (13)+3.B,(1 3)=C,(1 1 ) 




FIG.7B 



B R (9)-3.B«(13)-3.8 ( (13)gC R (11) 



-B,(9)-3.B R (13)O.B l (l3)=C,(11) 



B R (9)+3.B«(13)+3.B t (13)=C R (9) 



8,(9)-3.B a (13)+3.8,(13)=C,(9) | 
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A/B (0.1) (0,1) 



MP=0 
e1+e2 

e1-e2 

-(e3-e4) 

e3+e4 



MP=1 
PERM=0 

el*A.e3-B.e4 

e2+B.e3+A.e4\ 

e1-A.e3+B.e4 



MP=1 
PERM=1 

a e1-A.e3+B.e4 
-e2+B.e3+A.e4 
1 e1+A.e3-B.e4 



-e2+B.e3+A.e4' 1 e2+B.e3+A.e4 



FIG.11 
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MP=0 
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-(e3-e4) 
e3*e4 



MS=1 
MP=0 

<eUe2)+(e3+e4) 
e1-e2 
-(e3-e4) 
<e1+e2)-(e3+e4) 



MP-1 
PERM=0 



MP=1 
PERM=1 




MB (0.1) (0.1) 
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(x(0)+x(4)l+|x(2)+x(6)]=B R (0) 



x(0)-x(4)*A R (1)=B R (1) 



•tx(2)-x(6)]=.A R (3)=B,n) 



fx(0)+x(4)Hx(2)+x(6)j=B R (2) 



-«x(1 )*x(S)Hx(3)+x(7))1=-B B (6) I I 



x(1)-x(5)=A ft (5)=B R (S) 



-fx(3)-x<7)1=-A R (7)*B,(5) 



fx(1 )»x(5)l+(x(3)+x(7)l=B R (4) 




B R (0)^B B (4)=XA(0)=y(0) 



B R (1)+2.B R (5)+2.B,(5)=X R (1)=y(1) 



S,(1)-ZB»(5)+2.B,(5)=X.(1 )=y(2) 



B R (2)=X R (2)=y(3) 



■ B R (6)*X(2)*y(4) 



B B (1 )-2.B«(S)-2.B,(5)=X R (3)=y(5) 



-B»(1)-2.B R (5)+2.B,(5)=X l (3)=y(6) 



B R (0)-B ft (4)=X R (4)=y(7) 
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RAPPORT D'EXAMEN 
PRELIMINAIRE INTERNATIONAL 



Demande intemationale n° PCT/FR98/02636 



I. Base du rapport 

1 Ce rapport a ete redige sur la base des elements ci-apres (les feuilles de remplacement qui ont ete remises a 
I'office recepteur en reponse a une invitation faite conformement a I'article 14 sont considerees, dans ie present 
rapport, comme "initialement deposees" et ne sont pas jointes en annexe au rapport puisqu'elles ne contiennent 
pas de modifications.) : 

Description, pages: 

1 -38 version initiale 



Revendications, N°: 

1 ,3-42 version initiale 

2 re9ue(s) le 



04/01/2000 avec la lettre du 30/1 2/1 999 



Dessins, feuilles: 

1/17-17/17 



version initiale 



2. Les modifications ont entraine I'annulation : 

□ de la description, pages : 

□ des revendications, n os : 

□ des dessins, feuilles : 

3. □ Le present rapport a ete formule abstraction faite (de certaines) des modifications, qui ont ete considerees 

comme allant au-dela de I'expose de I'invention tel qu'il a ete depose, comme il est indique ci-apres 
(regie 70.2(c)) : 



4. Observations complementaires, le cas echeant : 
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V. Declaration motivee selon l article 35(2) quant a la nouveaute, lactivite inventive et la possibility 
d application industrielle; citations et explications a l appui de cette declaration 

1. Declaration 

Nouveaute Oui : Revendications 1-42 

Non : Revendications 

Activite inventive Oui: Revendications 1-42 

Non : Revendications 

Possibility d'application industrielle Oui: Revendications 1-42 

Non : Revendications 



2. Citations et explications 
voir feu i lie separee 



VIII. Observations relatives a la demande intemationale 

Les observations suivantes sont taites au sujet de la clarte des revendications, de la description et des 
et de la question de savoir si les revendications se fondent entierement sur la description : 

voir feuille separee 
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Concernant le point V 

Declaration motivee selon la regie 66.2(a)(ii) quant a la nouveaute, l activite 
inventive et la possibility d'application industrielle; citations et explications a 
I'appui de cette declaration 

1 ) II est fait reference aux documents suivants : 

D1:S. C. CHAN ET AL: "Software Implementation of 2D FFT on DSP 96002 
Digital Signal Processor" SIGNAL PROCESSING VI - THEORIES AND 
APPLICATIONS - PROCEEDINGS OF EUSIPCO-92 - SIXTH EUROPEAN 
SIGNAL PROCESSING CONFERENCE, 24 - 27 aout 1992, pages 1525-1528, 
XP000356534 BRUXELLES, BE 

D2 : US-A-5270953 

2) . L'invention concerne un procede de calcul de la transformed de Fourier rapide 

d'une suite de nombres reels ou d'une suite d'echantillons complexes conjugues. 

3) . Le probleme resolu par l'invention est la suppression des redondances des mises 

en memoire d' informations destinees a etre traitees par les circuits croisillons. 

4) . L'idee d'utiliser comme ordre de delivrance des echantillons de sortie I'inverse de 

I'ordre des echantillons de depart est connue dans le cas d'echantillons 
complexes (voir D1, page 1525, colonne de gauche, dernier paragraphe) est 
connue , de meme que de maniere generate les methodes de transformation "in- 
place" et la substitution en memoire des echantillons de sortie aux echantillons 
d'entree de rangs correspondants (D1, page 1526, colonne de gauche, 2e 
paragraphe). 

Par ailleurs, I'etat de la technique comprend les circuits de traitement separe des 
entrees pair/impaires et de calcul de somme des parties imaginaires (voir par 
exemple D2, revendication 3. 

5) . La solution consiste.dans le cas d'un signal numerique defini par une suite 

d'echantillons reels, a modifier les operations de calcul effectuees par les circuits 
croisillons aboutissant a la definition des echantillons de sortie et a substituer en 
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memoire les echantillons de sortie aux echantillons d'entree de meme rang 
correspondants de facon a delivrer les echantillons de sortie a la derniere etape 
de transformation dans I'ordre croissant de I'indice n, selon les relations : 
y(0) = Re[X(0)] 

y(n) = Re[X(tn+1 )/2] pour n impair et different de N-1 
y(n) = lm[X(n/2)] pour n pair et different de 0 
y(N-1) = Re[X(N/2)] 

ou les X(n) sont les echantillons complexes de la suite correspondant a la 
transformee de Fourier rapide. 

Dans le cas d'une suite d'echantillons complexes conjugues (revendication 2), la 
solution consiste a effectuer le calcul sur une suite de N echantillons reels de 
depart y(n) representative de la suite des N echantillons complexes X(n) 
conjugues deux a deux et definis par les relations ci-dessus. 

6). La solution n'est ni divulguee ni suggeree dans I'etat de la technique disponible. 



Concernant le point VIII 

Observations relatives a la demande intemationale 



7). Contrairement a ce qu'exige la regie 5.1 a) ii) PCT, la description n'indique pas 
I'etat de la technique anterieure pertinent expose dans les documents D1 et D2 et 
ne cite pas ces documents. 
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y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] pour n impair et 

different de N-l 

y(n) - Im[X(n/2)] pour n pair et different 

de 0 

y(N-l) = Re [X (N/2) ] 
oh les echantillons X(n), avec n € [0..N-1], designent 
les echantillons complexes de la suite correspondant a 
la transformee de Fourier rapide ou rapide inverse de 
la suite des echantillons de depart x(n). 

2. Proc§de de calcul de la transformee de Fourier 
rapide ou de la transformee de Fourier rapide inverse 
d'un signal numerique defini par une suite de N 
echantillons complexes X(n) conjugues deux a deux, 
caracteris6 en ce que le calcul est effectue sur une 
suite de N echantillons reels de depart y(n) 
representative de la suite d 1 echantillons complexes 
X (n) , avec N puissance de deux et n e [0..N-1], les 
echantillons de depart y(n) etant definis de la maniere 
suivante : 

y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+l) /2) ] pour n impair et 

different de N-l 

y(n) - Im[X(n/2)] pour n pair et different 

de 0 

y (N-l) « Re[X(N/2) ] 
en ce que le procede comporte des etapes successives de 
transformation pour transformer des echantillons 
d f entree en echantillons de sortie, une suite de N 
echantillons de sortie x(n) reels representatif s de 
cette transformee de Fourier rapide ou rapide inverse 
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§tant delivree a la derniere etape de transformation, 
l f ensemble des etapes de transformation etant effectu§ 
au moyen d'un unique jeu de circuits croisillons & 
plusieurs entrees et plusieurs sorties dont le mode 
operatoire est change s61ectivement a chaque etape de 
transformation, les echantillons d 1 entree et de sortie 
de chaque etape de transformation etant stock§s dans 
une memoire de stockage, 

et en ce que les echantillons de sortie d'un circuit 
croisillon se substituent dans la memoire de stockage 
aux echantillons d' entree de meme rang correspondants, 
de telle sorte que, si les echantillons de depart y(n) 
trait6s par la premiere etape de transformation sont 
classes dans I'ordre croissant de I'indice n, les 
echantillons de sortie x(n) sont delivres a la derniere 
etape de transformation dans I'ordre binaire inverse de 
I'indice n. 

3. Procede de calcul selon la revendication 1 ou 2, 
caracterise en ce que, a chaque etape de 
transformation, chaque circuit croisillon transforme 
des couples d' echantillon d ! entree, les rangs des 
echantillons d f entree d'un meme couple dans la suite 
des echantillons d' entree de ladite etape de 
transformation etant symetriques par rapport a un 
milieu entre les valeurs de rang extremes des 
echantillons d' entree transformers par ledit circuit 
croisillon . 

4. Procede de calcul selon la revendication 3, 
caracterise en ce qu'il comporte y-1 etapes de 
transformation E p avec p=log2 (N) et p e [0..p-2]. 



PEWLLE M0DIF1EE 
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deration pour apprecier la nouveaute et I'activite inventive de I'invention, objet de votre demande. 

Selon I'article R. 61 2-59 du code precite. vous disposez d'un delai de 3 mois a compter de la date de 
reception de ce rapport de recherche preliminaire pour y repondre par ecrit. Avant I'expiration de ce delai. 
ceiui-ci peut etre renouvele une fois sur votre requete. 



Suivant la categorie des documents cites, vous pouvez etre tenu a une obligation de reponse (par exemple. 
si le rapport de recherche preliminaire mentionne des documents de categorie X ou Y). Dans ce cas. un 
papilion rouge est appose sur cette lettre et le defaut de reponse entrainera le rejet de la demande. Dans 
le cas contraire, ce papilion est jaune. 

Dans tous les cas, il est de votre interet en elaborant votre reponse. de tenir compte de tous les 
documents cites. 



Selon les articles R.612-58 et R. 61 2-60 du code precite. votre reponse peut consister : 

- soit en de nouvelles revendications (en 3 exemplaires). Dans ce cas, vous devez signaler les 
changements apportes aux revendications initiales. Vous pouvez y joindre des observations qui mettent en 
evidence les caracteristiques techniques de ces nouvelles revendications qui echappent a I'opposabUite 
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1-3 



-3 
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1-42 

1-42 



DOMAINES TECHNIQUES 
RECHERCHES (lnt.CL.6) 



G06F 



Data d'achevement de la recherche 

25 aout 1998 



Examinateur 

Abram, R 



CATEGORIE OES DOCUMENTS CITES 

X : particulierement pertinent a lui seul 

Y : particulierement pertinent en combinaison avecun 

autre document de la me me categorte 
A : pertinent a rencontre d'au moins une revendication 

ou arriere-plan technologique general 
O : divulgation non-ecrite 
P : document intercalate 



T : theorie ou principe a ta base de Tinvention 

E : document de brevet beneficiant dune date anterieure 

a ta date de depdt et qui n'a ete public qu a cette date 

de dep6t ou qua une date posterieure. 
O : cite dans la demande 
L : cite pour d'au t res raisons 

&. : membre de la meme famille, document correspond ant 



•VYHSPAGEBUMKIWW 



ANNEXE AU RAPPORT DE RECHERCHE PRELIMINAIRE 

RELATIF A LA DEMANDE DE BREVET FRANQAIS NO. FA 553682 

FR 9715737 

La presente annexe indique les membres de la famille de brevets relatifs aux documents brevets cites dans le rapport de 
recherche prelim inaire vise ci-dessus. 

Lesdits membres sont contenus au fichier informatique de I'Office europeen des brevets a la date du 

Ues renseignements fournis sont donnes a titre indicatif et n'engagent pas la responsabilite de I'Orfice europeen des brevets. 

ni de 1'Administration franchise ^ qq 1998 



Document brevet cite 
au rapport de recherche 



Date de 
publication 



Membre(s) de la 
famille de breve t(s) 



Date de 
publication 



US 5270953 
US 4612626 



A 
A 



14-12-1993 
16-09-1986 



AUCUN 



AUCUN 



o 



Pour tout renseignement concernantcette annexe : voir Journal Orficiel de I'Orfice europeen des brevets, No.12/82 
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RAPPORT DE RECHERCHE INTERNATIONALE 
(article 18 et regies 43 et 44 du PCT) 



Reference du dossier du deposant ou 

du mandataire 

013770-PC. 


POUR SUITE voir ' a notification de transmission du rapport de recherche Internationale 
(formulaire PCT/IS A/220) et, le cas echeant, le point 5 ci-apres 

A DONNER 


Demands Internationale n° 
PCT/FR 98/02636 


Date du depdt international (jour/mois/annte) 

07/12/1998 


(Date de priorite (la plus ancienne) 
(jour/mois/ann£e) 

08/12/1997 


Deposant 

FRANCE TELECOM (SA) et al. 



Le present rapport de recherche international e, etabli par I'administration chargee de ta recherche Internationale, est transmis au 
deposant confornrtement a I'article 18. Une copie en est transmise au Bureau international. 

Ce rapport de recherche international comprend 3 feuilles. 

[X] II est aussi accompagne d'une copie de chaque document relatif a I'etat de la technique qui y est cite. 



1 . Base du rapport 

a. En ce qui concerne la langue, la recherche international a ete effective sur la base de la demande internationale dans la 
langue dans laquelle elle a ete deposee, sauf indication contraire donnee sous le m^me point. 



□ 



la recherche internationale a ete effectuee sur la base d'une traduction de la demande internationale remise a ('administration 



En ce qui concerne les sequences de nucleotides ou d'acides amines divulguees dans la demande internationale (le cas echeant) 
la recherche internationale a ete effective sur la base du listage des sequences : 
[ | contenu dans la demande internationale, sous forme ecrite. 

deposee avec la demande internationale, sous forme ddchiffrable par ordinateur. 

remis ulterieurement a ('administration, sous forme ecrite. 

remis ulterieurement a I'administration, sous forme dechiffrable par ordinateur. 



2. 
3. 



□ 
□ 
□ 
□ 

□ 

□ 
□ 



La declaration, selon laquelle le listage des sequences presente par ecrit et fourni ulterieurement ne vas pas au-dela de la 
divulgation faite dans la demande telle que deposee, a ete fournie. 

La declaration, selon laquelle les informations enregistrees sous forme dechiffrable par ordinateur sont identiques a cedes 
du listage des sequences presente par ecrit, a ete fournie. 

II a ete estime que certalnes revendications ne pouvaient pas faire 1'objet d'une recherche (voir le cadre I). 
II y a absence d'unite de ('Invention (voir le cadre II). 



4. En ce qui concerne le litre, 

[X"| le texte est approuve tel qu'il a ete remis par le deposant. 

| | Le texte a ete etabli par I'administration et a la teneur suivante: 



5. En ce qui concerne Pabrege, 

le texte est approuve tel qu'il a ete remis par le deposant 

□ le texte (reproduit dans le cadre Ml) a ete etabli par I'administration conformement a la regie 38.2b). Le deposant peut 
presenter des observations a I'administration dans un deiai d'un mois a compter de la date d'expedition du present rapport 
de recherche internationale. 

6. La figure des dessins a publier avec I'abrege est la Figure n° 8 

[X] suggeree par le deposant. | | Aucune des figures 

| | parce que le deposant n'a pas suggere de figure. n eSt ^ P u * )l ' er - 

| | parce que cette figure caracterise mieux ('invention. 
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A. CLASSEMENT DE L'OBJET DE LA DEMANDE 

CIB 6 G06F17/14 



Solon ta classification Internationale des brevets (CIB) ou a la fois selon la classification natron ale et la CIB 



B. DO MAIN ES SUR LESQUELS LA RECHERCHE A PORTE 



Documentation minimale consuttee (systeme de classification sutvl des symboles de classement) 

CIB 6 G06F 



Documentation consulted autre que la documentation minimale dans la mesure ou ces documents relevent des domaines sur lesquels a porte la recherche 



Base de donnees electronique consultee au cours de la recherche Internationale (nom de la base de donnees, et si realisable, termes de recherche utilises) 



C. DOCUMENTS CONSIDERES COMME PERTINENTS 



Categorie ° Identification des documents cites, avec, le cas echeant, IMndication des passages pertinents 



no. des revendi cat ions visees 



S. C. CHAN ET AL: "Software Implemental on 
of 2D FFT on DSP 96002 Digital Signal 
Processor" 

SIGNAL PROCESSING VI - THEORIES AND 
APPLICATIONS - PROCEEDINGS OF EUSIPCO-92 - 
SIXTH EUROPEAN SIGNAL PROCESSING 
CONFERENCE, 24 - 27 aoGt 1992, pages 
1525-1528, XP000356534 
BRUXELLES, BE 

voir le document en entier 

S. L. JOHNSON: "Cooley-Tukey FFT on the 

Connection Machine" 

PARALLEL COMPUTING, 

vol. 18, no. 11, novembre 1992, pages 

1201-1221, XP000330876 

AMSTERDAM, NL 

voir le document en entier 

-/-- 



1-3 



4-42 
1-3 



4-42 



Voir la suite du cadre C pour la fin de ta liste des documents 



Les documents de families de brevets sont indiques en annexe 



0 Categories spectales de documents cites: 

"A" document definissant Tetat general de la technique, non 
considere com me particulierement pertinent 

"E" document anterieur. mais publie a la date de depdt international 
ou apres cette date 

"L" document pouvant jeter un doute sur une revendication de 
priorite ou cite pour determiner la date de publication o"une 
autre citation ou pour une raison speciale (telle qu'indiquee) 

"O" document se referant a une divulgation orale, a un usage, a 

une exposition ou tous autres moyens 
"P" document publie avant la date de depdt international, ma is 

posterieurement a la date de priorite revendiquee 



T" document uft6rieur publie apres la date de depot international ou la 
date de priorite et n'appartenenant pas a I'etat de la 
technique pertinent, mais cite pour comprendre le principe 
ou la theorie constituant la base de I'invention 

"X" document particulierement pertinent; ttnven tion revendiquee ne peut 
etre consideree comme nouvelle ou comma impliquant une activite 
inventive par rapport au document considere isolement 

"Y" document particulierement pertinent; finven tion revendiquee 

ne peut etre consideree comme impliquant une activite inventive 
lorsque le document est associe a un ou ptusieurs autres 
documents de memo nature, cette combinaison etant evidente 
pour une personne du metier 

document qui fart partie de la meme famille de brevets 



Date a laquelle la recherche Internationale a ete eff ectivement achevee 



30 avril 1999 



Date d' expedition du present rapport de recherche intemationate 



11/05/1999 



Nom et adresse postale de I' administration chargee de la recherche Internationale 
Office European des Brevets, P.B. 5818 Patentlaan 2 
NL - 2280 HV Rtjswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Fonctionnaire autorise 



Abram, R 
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Categorte 0 Identification des documents cites, avec.le cas echeant, Tlndlcatlondes passages pertinents 



no. das rev© indications vis^es 



US 5 270 953 A (S. A. WHITE) 

14 decembre 1993 

voir le document en entier 
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16 septembre 1986 

voir le document en entier 
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Document brevet cite* 


Date de 


Membre(s) de la 


Date de 


au rapport cle recherche 


publication 


famitle de brevet(s) 


publication 



US 5270953 



14-12-1993 



AUCUN 



US 4612626 A 16-09-1986 AUCUN 
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PATENT COOPERATION TRE^TV 

PCT 

INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



(PCT Article 36 and Rule 70) 


Applicant's or agent's file reference 
013770-PC. 


r*o itt dtufd rTIrtw See Notification of Transmittal of International 
FOR FURTHER ACTION PreIirninary Examination Report (Form PCT/IPE A/4 1 6) 


International application No. 


International filing date (day/month/year) 


Priority date {day/month/year) 


PCT/FR98/02636 


07 December 1998 (07.12.98) 


08 December 1997 (08.12.97) 


International Patent Classification (IPC) or national classification and IPC 




G06F 17/14 






Applicant 


FRANCE TELECOM (S.A.) 





This international preliminary examination report has been prepared by this International Preliminary Examining 
Authority and is transmitted to the applicant according to Article 36. 



This REPORT consists of a total of 



. sheets, including this cover sheet. 



<k7\ This report is also accompanied by ANNEXES, i.e., sheets of the description, claims and/or drawings which have 
^ been amended and are the basis for this report and/or sheets containing rectifications made before this Authority 
(see Rule 70.16 and Section 607 of the Administrative Instructions under the PCT). 



These annexes consist of a total of 



sheets. 



This report contains indications relating to the following items: 
Basis of the report 
Priority 

Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 
Lack of unity of invention 

Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

Certain documents cited 

Certain defects in the international application 

Certain observations on the international application 



1 


I2SI 


II 


□ 


III 


□ 


IV 


□ 


V 


[El 


VI 


□ 


VII 


□ 


VIII 





Date of submission of the demand 

05 July 1999 (05.07.99) 


Date of completion of this report 

17 February 2000 (17.02.2000) 


Name and mailing address of the IPEA/EP 
Facsimile No. 


Authorized officer 
Telephone No. 



Form PCT/IPEA/409 (cover sheet) (January 1994) 



n 



] 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



International application No. 

PCT/FR98/02636 



I. Basis of the report 



1 . This report has been drawn on the basis of (Replacement sheets which have been furnished to the receiving Office in response to an invitation 
under Article 14 are referred to in this report as ' originally filed" and are not annexed to the report since they do not contain amendments.): 



[ [ the international application as originally filed. 

J^Xj the description, pages *"38 

pages 

pages 

pages 



as originally filed, 
filed with the demand, 
, filed with the letter of 
, filed with the letter of 



the claims, 



Nos. 
Nos. 
Nos. 
Nos. 
Nos. 



1, 3-42 



as originally filed, 

as amended under Article 19, 

filed with the demand, 
filed with the letter of 
, filed with the letter of 



30 December 1999 (30.12.1999) 



the drawings, 



1/17-171/17 



sheets/fig 

sheets/fig 

sheets/fig 

sheets/fig 

2. The amendments have resulted in the cancellation of: 

1 1 the description, pages 

1 I the claims, Nos. 

i I the drawings, sheets/fig 



, as originally filed, 
, filed with the demand, 
, filed with the letter of 
, filed with the letter of 



2 I I This report has been established as if (some of) the amendments had not been made, since they have been considered 
I — I t0 g 0 b e y 0nc i the disclosure as filed, as indicated in the Supplemental Box (Rule 70.2(c)). 



4. Additional observations, if necessary: 



Form PCT/IPEA/409 (Box I) (January 1994) 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



^^fiational application No. 
PCT/FR 98/02636 



V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 



Statement 
Novelty (N) 

Inventive step (IS) 
Industrial applicability (IA) 



Claims 
Claims 

Claims 
Claims 

Claims 
Claims 



1-42 



1-42 



1-42 



YES 
NO 
YES 
NO 

YES 
NO 



Citations and explanations 

1) Reference is made to the following documents: 

Dl : S.C. CHAN ET AL : "Software implementation of 2D 
FFT on DSP 96002 Digital Signal Processor" SIGNAL 
PROCESSING VI - THEORIES AND APPLICATIONS - 
PROCEEDINGS OF EUSIPCO-92 - SIXTH EUROPEAN SIGNAL 
PROCESSING CONFERENCE , 24-27 August 1992, pages 
1525-1528, XP000356534 Brussels, BE 



D2: US-A-5270953 



2) The invention relates to a method for computing the 
fast Fourier transform of a series of real numbers 
or a series of complex conjugated samples. 



3) The problem solved by the invention is eliminating 
the redundancy of storing data to be processed by 
crossbar circuits . 



4) The idea of using as the output sample delivery 

order, the reverse of the initial sample order, is 
known in the case of complex samples (see Dl, page 
1525, left-hand column, final paragraph) and, 
generally speaking, the in-place methods of 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



^j^iational application No. 
PCT/FR 98/02636 



transformation and the substitution in memory of 
output samples for input samples having 
corresponding ranks are also known (Dl, page 1526, 
left-hand column, second paragraph) . 
Moreover, the prior art includes circuits for 
processing the odd/even inputs separately and 
computing the sum of the imaginary portions (see, 
for example, D2, Claim 3) . 



5) In the case of a digital signal defined by a series 
of real samples, the solution consists of modifying 
the computing operations carried out by the crossbar 
circuits leading to the definition of output 
samples, and substituting in the memory the output 
samples for the corresponding input samples having 
the same rank so as to deliver the output samples to 
the final stage of transformation in the ascending 
order of index n, according to the relationships: 

y (0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2] when n is odd and different 

from N-l 

y(n) = Im[X(n/2)] when n is even and different 

from 0 

y (N-l) = Re[X(N/2) ] 

wherein the X(n) are complex samples from the series 
corresponding to the fast Fourier transform. 
In the case of a series of complex conjugated 
samples (Claim 2), the solution consists of carrying 
out a computation on a series of N real initial 
samples y(n) representative of the series of N 
complex samples X(n) conjugated in pairs and defined 
by the above relationships. 

6) The solution is not disclosed in, or suggested by, 
the available prior art. 
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onal application No. 
PCT/FR 98/02636 



VIIL Certain observations on the international application 



The following observations on the clarity of the claims, description, and drawings or on the question whether the claims are fully 
supported by the description, are made: 

7) Contrary to the requirements of PCT Rule 5.1 

(a) (ii), the description does not indicate the 
relevant prior art disclosed in documents Dl and D2, 
nor does it cite this document. 
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